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La dilatation thermique et la chaleur spécifique à volume constant. Loi de Grüneisen. 


Résumé. — Selon Grüneisen, le coefficient de dilatation et la chaleur spécifique à volume cons- 
tant resteraient dans le même rapport à toute température. En évaluant la dilatation thermique 
à partir de l’énergie potentielle du milieu cristallin, on aboutit à la conclusion que la loi de Grü- 
neisen n’est qu’approximative ; à peu près en accord avec les faits aux températures moyennes (tant 
que la chaleur spécifique à volume constant garde sensiblement la même grandeur), mais nota- 
blement en défaut aux températures basses et élevées. 

D’autre part, on trouve que la constante de Grüneisen est déterminée principalement par le 
rapport entre les grandeurs absolues des dérivées troisièmes et des dérivées secondes de l’énergie 
potentielle du milieu cristallin, et non par la variation de la fréquence des oscillations thermiques 
en fonction de la dilatation. Enfin, on établit que, si le milieu cristallin est libre de toute contrainte, 
les dérivées secondes de son énergie potentielle, les fréquences des oscillations thermiques et les 
coefficients d’élasticité sont, en gros, des fonctions exponentielles décroissantes de la dilatation 
thermique. Pris 


Abstract. — According to Grüneisen, the expansion coefficient and the specific heat at constant 
volume would remain in the same ratio at any temperature. The estimation of the thermal ex- 
passion from the potential energy of the crystalline medium leads to the conclusion that Grüneisen’s 
law is only approximate ; it is nearly in agreement with facts at medium temperatures (as long 
as the specific heat at constant volume remains about thesame magnitude), but at low and high 
temperatures it is greatly at variance. 

Moreover, we find that the Grüneisen constant is mainly determined by the ratio between the 
absolute magnitudes of the third and second derivatives of the potential energy of the crystalline 
medium, and not by the variation of the frequency of the thermal oscillations in relation to the 
expansion. Lastly we show that, if the crystalline medium is free from constraint, the se- 
cond derivatives of its potential energy, the frequencies of the thermal oscillations and the elastic 
coefficients are roughly decreasing exponential functions of the thermal expansion. 
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A. Le coefficient de dilatation et la chaleur spéeci- 
fique à volume constant. Loi de Grüneisen (1). — 
Supposons les tensions secondaires négligeables et 
les facteurs Es. (24) invariables, tous égaux, de 
grandeur £ (41-42). En ce cas, la tension ther- 
mique 7 (12) est proportionnelle à l’énergie totale U 
des oscillations fondamentales : 


EU 
N uv 


(:) Les formules numérotées de 1 à 31, les références 
de [1] à [5],les figures 1 et2 se trouvent dans la première 
partie (p.451) ; les formules numérotées de 32 à 53, les 
références de [6] à [10], les figures 3, 4, 5, 6, 7, dans la 
seconde partie (p. 465). 


1 : U = À Wgr. 
ST 


Et le coefficient de dilatation linéaire devient 


Ek | 
4 je on à (54) 
H=x2 [re (Si) d Log de 


Admettons aussi que le coefficient € (8) est inva- 
riable, et que 34 peut être tenu pour nul par rap- 
port à uC. Il vient : 


(e 


Ek Cov 
Ne dm = UU 


Co eSt la chaleur spécifique à volume constant qui 
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est prise par les oscillations fondamentales ; et 
selon notre hypothèse, & Se confond avec la cha- 
leur spécifique définie par la théorie classique : 


cg = (kIN) 2 B$r. 


Si toutes les conditions que nous venons de poser 
étaient satisfaites, le coefficient de dilatation li- 
néaire serait proportionnel à 
fique co qui revient aux oscillations fondamen- 


tales : 


(55) 


afco = EC. 


Mais les facteurs Ë$. (24) sont inégaux.et varient 
légèrement avec la température ; le coefficient € di- 
minue quand le milieu cristallin se dilate ; aux tem- 
pératures élevées les tensions secondaires 7, (28), 
T3 (29), ... contribuent à la dilatation ; et la 
tension 7, est du second degré, la tension 7, du 
troisième degré par rapport aux énergies Ws. (3) 
des oscillations fondamentales. Donc le rap- 
port «fc Varie quand la température s’élève ; et sa 
variation est considérable aux basses et aux hautes 
températures. Car si la température est basse, 
quand elle s’élève les fonctions 65, (4) ne restent pas 
entre elles dans des rapports constants, et quand la 
température est haute, les tensions secondaires 
prennent part à la dilatation. 

Toutefois, aux températures moyennes — qui se 
rapportent à la branche BC de la courbe OD 
(fig. 7, J. Physique Rad. Août-Sept 1961, 22. 
fase. 8-9, 470). — toutes les fonctions 64, (4) sont 
sensiblement égales, car toutes atteignent sensible- 
ment l’unité ; et les tensions secondaires restent 
encore négligeables par rapport à la tension princi- 
pale ; en conséquence, l’égalité (56) est à peu près 
satisfaite. D’où l’on déduit, compte tenu de la 
relation (8) entre le coefficient € et le module de 
compressibilité isotherme x, une nouvelle égalité : 


TRUSSS 


(56) 


Re. (57) 
où V est le volume spécifique : 
1 
Pole (B = 30, =5) ICS) 


Le facteur moyen Ë (25) reste à peu près constant 
quand le milieu cristallin se dilate. Donc, aux tem- 


B V 
pératures moyennes, le rapport 1. est aussi 
00 


approximativement constant ; il prend les gran- 
deurs suivantes (déduites des facteurs moyens Ë 
évalués précédemment (42)). 
Aj Fecx 
OT er 

mn ? ) 


Cu 
4,9 


B. La théorie et la loi de Grüneisen. —— Mais la 
chaleur spécifique, «, à volume constant, ne se 
réduit pas à la chaleur spécifique co (55) qui est 
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prise par les oscillations fondamentales. Elle com- 


prend en plus trois termes : l’un c,, d’origine élec- 


tronique, l’autre c», d’origine magnétique, impor- - 


tant seulement si le cristal est ferromagnétique, le 
troisième c,, revient aux oscillations secondaires : 


Co = Cov + Cev + Cmv + Csv. 


Grüneisen a mesuré sur les cristaux cubiques 


envisagés, à des températures moyennes, les coef- ” 


ficients de dilatation et les chaleurs spécifiques [14]. 


Il trouve que le rapport (B/c,) est sensiblement … 
constant. C’est qu'aux températures moyennes, 


toutes les fonctions 64, (4) sont à peine moindres 
que l’unité ; donc à ces températures 


CR a 
Cy He [1 + (Cev + Cmo + Csv) [Cor] 


Et, quand la température s’élève, le rapport . 


(Co + Cmo + Cæ) Co Augmente, et cette augmen- 
tation compense, au moins en partie, la diminution 


du coefficient C (8). Mais aux températures basses … 


et élevées, les chaleurs spécifiques Ce, Cms, Ne Crois- 
sent pas avec la température suivant la même loi 
que le coefficient de dilatation $ (1). Quant à la 
chaleur spécifique c,, négligeable aux basses et aux 


moyennes températures, elle ne reste pas, non plus, 


dans un rapport constant avec le coefficient de dila- 
tation 6 lorsque la température du cristal devient 
proche du point de fusion. Donc aux températures 
basses et élevées, le rapport B/c, varie notablement. 

Pour rendre compte de ses expériences Grüneisen 
se réfère à la thermodynamique [11] et fait deux 
hypothèses : 

19 Les composantes harmoniques de l’agitation 
thermique sont indépendantes les unes des autres ; 


done chacune a la même énergie libre qu’un oscil- « 


lateur harmonique de même fréquence. Ainsi l’éner- 
ge libre de l’agitation thermique est égale à la 


somme, S, des énergies libres prises par les oscil-. 


lations harmoniques qui forment cette agitation (2) 


ia __ ho ño 
s—> [ke Log. (1 ES re) +]. (60) 
Et l’énergie libre du cristal s’élève à 
A=W+S (61) 


VW, (1) étant l’énergie potentielle du cristal sans 
agitation thermique, c’est-à-dire : les atomes au 
repos dans leurs positions moyennes (3). 


(*) D’après les théories présentement admises et d’après 
la calorimétrie, la chaleur spécifique ce est grossièrement 
proportionnelle à la température absolue, et la chaleur 
spécifique Cm» à peu près en raison directe de la dérivée 
do? fd®, 6 étant l’aimantation spécifique. 


(*) Je définis ici (formule 60) l’énergie libre de l’agitation 


thermique par une expression explicite ; Grüneisen rap- 
pare cette énergie à la température caractéristique de 
ebye. 
(5) Formule (10), page 17-18, Handbuch der Physik, 
Band X. 


(59) . 


| Nos 


| 29 Quand le volume spécifique, V(58), du milieu 

cristallin augmente, donc quand le cristal se dilate, 
les fréquences vs des oscillations harmoniques 
faites par les atomes s’abaissent, et les dérivées 
dLog.og 
d Log.v 
température. Subsidiairement, Grüneisen suppose 
que toutes les dérivées d Log.o,./d Log.v scnt 
d Log. 


sont constantes, donc indépendantes dela 


sensiblement égales, de grandeur RTE ( 
Et de ces données, il déduit l’égalité : 
, BV  dLog.o 
Cox dLog.v (a) 


Toutefois, dans cette égalité, Grüneisen fait état 
non de la chaleur spécifique globale c, (59), mais de 
la chaleur spécifique classique co, (b3), celle qui est 
prise par les oscillations harmoniques des atomes 
supposées indépendantes (?). Donc le premier 
membre de l’égalité (57) se confond avec le premier 
membre de l’égalité (62) obtenue par Grüneisen. 


Mais je trouve que le rapport 1e 5 ê est égal à 5 et 
Co X 2 
AND UOP. Ge 
non au module de la dérivée Tir (PL 


Les composantes harmoniques de l'agitation 
thermique ne sontindépendantes qu’autant qu’elles 
développent entre les atomes des forces de rappel 
rigoureusement conformes à la loi de Hooke. Cela 
implique que l’énergie potentielle du cristal n’a pas 
de dérivées de rang supérieur à deux, c’est-à-dire 
que ses dérivées troisièmes et suivantes sont toutes 


(?) En première étude [11] Grüneisen admet, comme 
Einstein, que l’agitation thermique est formée d’oscillations 
harmoniques indépendantes qui ont toutes la même fré- 
quence. En outre, il suppose les atomes du milieu cristallin 
unis par des forces centrales, et rapporte la dérivée 
d Log.os-/d Log.v à de telles forces de cohésion. 

(?) C’est qu’au temps où Grüneisen développa sa théorie, 
les chaleurs spécifiques d’origine électronique et magné- 
tique étaient inconnues, et l’on confondait alors la chaleur 
spécifique co (55) avec la chaleur spécifique globale à 
volume constant € (59). 

(5) Précisément Grüneisen fait d’abord état de déri- 
vées partielles (d Log.ws-/àd Log.v)e ; puis il admet impli- 
citement que ces dérivées partielles expriment la variation 
que les pulsations ws- subissent lorsque le cristal se dilate 
libre de toute contrainte. Cela est contraire aux faits. La 
dilatation thermique est entraînée par une élévation de la 
température ; donc les dérivées des pulsations ws- rela- 
tives à la dilatation thermique se rapportent à la fois à 
une augmentation de volume et à une élévation de la 
température, et sont par là des dérivées totales. 

Surtout les tensions thermiques ne peuvent être rappor- 

ù Log. 28) ; 
ù Log.v /o 
Car, si le milieu cristallin est maintenu à volume constant, 
nulle variation des pulsations os. ne peut être produite 
par une variation de volume, il n’existe donc plus de 
dérivées (d Log.os-/d Log.v)e, tandis que les tensions 
thermiques subsistent à volume constant et sont à tempé- 
rature égale plus fortes que dans le milieu cristallin libre 
de toute contrainte. 


tées, comme le fait Grüneisen, aux dérivées 
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nulles. Si l’énergie potentielle est telle, ses dérivées 
secondes C% (2) sont constantes, et les pulsa- 
tions w$, sont constantes aussi (21) : elles ne varient 
point avec le volume spécifique Y (58) du milieu 
cristallin, et les dérivées d Log.os./d Log.v sont 
donc nulles ; en outre, une élévation de la tempé- 
rature n’entraîne nulle dilatation du cristal. 

Au contraire, si les pulsations w,, varient avec le 
volume spécifique V du milieu cristallin, -l’énergie 
potentielle du cristal admet des dérivées de rang 
supérieur à deux (4). Et les forces de rappel entre 
les atomes n’obéissent plus à la loi de Hooke. En 
conséquence, l’agitation thermique n’est plus la 
résultante d’oscillations harmoniques indépen- 
dantes ; donc, l’énergie libre de cette agitation n’est 
plus conforme à l’expression classique (60), et 
l’énergie libre du cristal ne répond plus à l’expres- 
sion (61). 

Grüneisen fait donc deux hypothèses incompa- 
tibles : l’une excluant l’autre. La relation qu’il 
obtient (62) ne peut être admise en raison de cette 
contradiction. 


C. Les fréquences des oscillations fondamentales 
et la dilatation thermique. — A température crois- 
sante, la structure électronique du milieu cristallin 
évolue. Et l’évolution est multiple. 

Supposons encore le milieu cristallin libre de 
toute contrainte. Lorsque la dilatation grandit, la . 
densité électronique diminue, et les énergies asso- 
ciées aux états des électrons varient. Cependant la 
diminution de la densité électronique est infime 
dans le corps central des atomes, elle ne prend une 
grandeur notable qu’à leurs confins ; de même, 
parmi les énergies associées aux états des électrons, 
seules varient sensiblement celles qui sont les plus 
faibles en grandeur absolue — celles qui concernent 
les états pris par les électrons de valence ou de 
conductibilité —. 

D’autre part, quand l’agitation thermique s’am- 
plifie, les échanges d’énergie se multiplient entre 
les ondes élastiques et les électrons faiblement liés. 
En conséquence, des électrons faiblement liés, en 
nombre croissant avec la température, sont portés 
à des niveaux d’énergies plus élevées que les niveaux 
où ils se trouvaient primitivement. En outre, si le 
cristal est ferromagnétique, l’agitation thermique 
des moments magnétiques portés par les atomes 
croît avec la température (5). 

La structure électronique évolue même s1 le 
milieu cristallin est maintenu à volume constant. 
En ce cas, la densité électronique et les énergies 
prises par les électrons ne varient que dans une 


(t) L'énergie potentielle du milieu cristallin a été le sujet 
de nombreuses études: Toutes les fonctions proposées pour 
la définir sont indéfiniment dérivables, en particulier celle 
que Grüneisen prend en compte. 

(t) Je ne fais pas état des oscillations de plasma ; car, 
selon les conceptions admises, elles gardent sensiblement la 
même énergie à toute température. 
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proportion minime,même si l’élévation de la tempé- 
rature est ample, de quelques centaines de degrés 
centésimaux ; mais le transfert d’électrons faible- 
ment liés vers des niveaux d’énergies plus élevées 
n’est sans doute pas fortement diminué. 

Cette évolution de la structure électronique, à 
température croissante, réduit, en grandeur absolue, 
l'énergie potentielle fondamentale W,(O) (1) du 
cristal, et, dans l’ensemble, toutes les dérivées de 
cette énergie : secondes C%, troisièmes DE, etc. (?). 
Elle a donc pour effet d’abaisser les fréquences des 
oscillations fondamentales ; et cet abaissement de 
fréquence se produit encore, mais devient minime 
lorsque le volume du cristal reste constant. 

Restreignons encore notre étude aux cristaux 
formés par les métaux, et envisageons de nouveau 
le milieu cristallin abandonné à lui-même. 

Si à la température O, les dérivées secondes, troi- 
sièmes, quatrièmes, de l’énergie potentielle 
%,(O) (1-2), sont respectivement C#(O), DÉ"(O), 
ESY(O) ..., à la température O@ + dO, la dérivée 
seconde C% devient (?) : 


(20) 
CÉ(E + 40) = CHE) + _ d® 


—Y DÉY(O) b, « d@ 
BY 


LUE [O) à, csla dO)? — -.. 


- 63 
2 (63) 


d 
D’où l’on déduit, les indices © omis : 


d Log.Cg&f 1 Fe Û ee (64) 
v ù Log.v /®@ 


dLog.v BC \30 
et 
Si Er VERS 1 By 
d Log.v Js LaÈ Dai br. (65 
À quelques rares expressions près, les deux déri- 
vées 
Les /? Log. Ca 
dO /» Ü 3 Log.v /e 


sont négatives (2). La première (dC%/d0),, prise à 


1) J'ai défini précédemment [5] les dérivées secondes, 
Ca” (®), qui se rapportent au cristal dilaté, en fonction des 
dérivées CS, DÉÉY, au zéro absolu de la température. Pour 
obtenir la relation (63) il suffit de prendre dans l'expression 
trouvée (10, 15-I) la température initiale égale à © et la 
dilatation égale à « d@ et d’ajouter la dérivée (dC% 20), 
que j'avais négligée. 

(*) Si l’on retient seulement les énergies potentielles bi- 
atomiques on a:b—a—m—bphpet 


2 Log. CE? 1 
(ee © 32 Dita, 
rl ue les cristaux cubiques formés par les métaux les 
rivées : 
d Log. 
d Log.v 
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volume constant par rapport à la température, est 
déterminée par le transfert des électrons de conduc-. 
tibilité vers des niveaux d’énergies élevées, -et, de 
plus, si le cristal est ferromagnétique, par l’ampli- 
fication apportée à l’agitation thermique des mo- 
ments magnétiques ; elle relève donc des valeurs 


spécifiques électronique et magnétique. La seconde … | 


(ic) 
dérivée, (ET) , prise à température cons- 
d Log.v }o 

tante par rapport à la variation de volume, est … 
fonction essentiellement de la différence, c, — & … 
entre la chaleur spécifique sous pression constante, 
Ch et la chaleur spécifique à volume constant, «. 

À température croissante, les dérivées troisièmes 
D#Y (2) diminuent, dans l’ensemble, en grandeur 
absolue, plus vite que la dérivée seconde C% (2). 


: Par là, le numérateur qui figure au second membre … 


de l’expression (65) décroît plus rapidement que le 
dénominateur de la même expression. Mais cet 
excès de décroissance est annulée approximati- 
vement par l’amplification des transiations 
b — m— q. Comme chaque facteur Es (24), et 
pour une raison semblable, toute dérivée 


É Lee.) 
ù Log.v /e 
reste à peu près constante du zéro absolu de la 
température au point de fusion (4). 

D'autre part, quand la température s’élève, les 
transferts des électrons faiblement liés se multi- 


(20) 
a 


dC ) Si 
0 /, #4 
mentent dans l’ensemble, tandis que les modules des 
produits Bo, C& diminuent presque tous (f) ; done 
les modules des dérivées 
1 = 
BC | 30 /, 


augmentent en très forte majorité. 

Sauf au voisinage du zéro absolu, pour une même 
élévation de la température, l’énergie potentielle 
VW,(0) (1-2) des cristaux envisagés est davantage 


plient, et les modules des dérivées 


sont probablement toutes négatives. Mais il n’est pas 
exclu que, pour certains cristaux, surtout si leur degré de 


d Log. Ci k 
symétrie est faible, quelques dérivées + soient 
positives. d Log.v 


Ho: CES 

4) La d û © 
(6) È érivée ( es 1h 
tante si les forces de cohésion entre les atomes du milieu 
cristallin étaient strictement conformes à la loi de Coulomb ; 
c’est-à-dire si entre deux atomes séparés par une distance r 
l’énergie de cohésion était cfr, e étant une constante. La 
dérivée (RE Log . Cé° 

à Log.v 
de cohésion était c/r?,. 

(5) Pour une même élévation de la température la dimi- 
nution relative de [Ca®] est, en général, beaucoup plus 
forte que l’augmentation relative du coefficient B(O), sauf 
aux basses températures. 


serait rigoureusement cons- 


je resterait encore constante si l’énergie 
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affaiblie, en grandeur absolue, par la dilatation que 
par le transfert des électrons de conductibilité vers 
des niveaux d’énergies plus élevées. D'autre part, la 
dilatation en écartant les atomes, réduit direc- 
tement et amplement les modules des dérivées 
secondes C%, Au contraire, le transfert des électrons 
faiblement liés qui diminue notablement l’énergie 
potentielle n’affaiblit que fort peu le module de ses 
dérivées. Les basses températures écartées, il est 
dCS : 
scæ (50 FORT, 
dans l’ensemble, petites en grandeur absolue par 
d pe) 
d Log.v Jo | 
Toutefois, les cristaux cubiques formés par les 
métaux de transition (M, Fu, N;, ...) ont un 
comportement particulier. Leur chaleur spécifique 
d’origine électronique est forte, probablement plus 
grande à toute température que la différence 
€, — entre la valeur spécifique sous pression cons- 
tante et la chaleur spécifique à volume constant. 
Donc chez ces cristaux, les modules des dérivées 
1 ee 
BCS | 20 
doute, du même ordre de grandeur que les modules 
ter É rec) 

des dérivées | 21) . 

ùd Log.v }e 
Tout compte fait, dans leur ensemble, c’est-à- 
dire en très forte majorité, les dérivées totales 
d Log.C£® 
d Log.v 
la température, mais les accroissements relatifs 
(ou les décroissances relatives) de leurs modules 


donc à prévoir que les dérivées 


rapport aux dérivées ( 


) , tout en restant plus petits, sont, sans 
v 


(64) croissent, en grandeur absolue, avec 


restent faibles entre le zéro absolu de la tempéra- 


ture et le point de fusion, de l’ordre de quelques 


- centièmes, même pour les cristaux formés par les 


métaux de transition. 
On trouve aussi : 
dLog.Dr 1 (oDS" 
d'Log-r. gDafr hO) 


1 + E°8rs 
cô 


= 3D* abc CG. (66) 


Et, du zéro absolu de la température au point de 
fusion, les dérivées 


d Log.DEfT [a Log.v 
comme les dérivées 
à Log.C£/a Log.v 


et pour des causes semblables, restent à peu près 
constantes. 
D'où découlent les expressions (18) où 


en ( ses Tee Sr 
sie ADS 0 /e _dLog.v /e 
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: (FE): nÊTS 1 (5) 
v v 


aps SN fe 
# pce (20 pDEf | 36 


Les oscillations fondamentales sont entretenues 
par des forces de rappel qui, à des écarts minimes 
près, restent linéaires par rapport aux dérivées 
secondes C8 (1-2). Négligeons ces écarts, et suppo- 
sons les vecteurs c$, constants, condition qui est 
aussi sensiblement réalisée. Compte tenu de la rela- 
tion (21), nous obtenons : 


d Log. os 
d Log.v 


2 Causr case D CP (SP + 2788 d) 
(Pa + na 0)8 
>> ce (Ps + n%8)à 


a 


sin? x.Sa 


»> Ca, Sr CR,ST sin? x Sa 
af 

(67) 
no cut: 


Quand la température s’élève, les modules des 
dérivées secondes CS ne décroissent pas unifor- 
mément ; donc chaque dérivée 


d Log.ws-[d Log.v 


varie, et suivant une loi particulière (déterminée 
par S et +). Mais seules importent les dérivées 
secondes C% qui portent sur des atomes voisins, et 
les exposants correspondants (o% + n%9)5 n’attei- 
gnent que quelques dixièmes au point de fusion. 
De ce fait, entre le zéro absolu de la température 
et le point de fusion, les dérivées 


d Log.ws-]d Log.v 


ne varient, en grandeur absolue, que dans une 
faible proportion (de quelques centièmes par rap- 
port à leur module au zéro absolu de la tempéra- 
ture) ; et la dérivée moyenne (négative) 


d Log. 1 fs d Log.os 


d Log.v SN d Log.v (68) 


varie encore moins. 

En conséquence, le facteur E/u (1) (57) et le 
module de la dérivée d Log. o/d Log.v (62) restent, 
à température variable, dans un rapport à peu près 
constant. Mais, comme les modules des dérivées 
secondes C et troisièmes DST ne diminuent pas 
uniformément à température croissante, il n’existe 


pas de relation simple entre le facteur et la déri- 


vée d Log.w/d Log.v. 
De même que la dérivée moyenne MER 
{ Ne d Log. ? 
(1) Le facteur À se confond avec la constante y de Grü- 
neisen [11]. 
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et pour la même raison, du zéro absoïu à la tempé- Toutes les dérivées (64-66-67-69) dont je viens 
rature de fusion, les dérivées (négatives) de faire état se rapportent au milieu cristallin libre. 
de toute contrainte, soit à l’évolution de ce milieu 


Dhs, Et re cree représentée par la courbe ODL (fig. 8) qui décrit la 
ee En nl See dilatation cubique 5(@) (16) en fonction de la tem- 
d Log.r (69) pérature, sous une pression nulle. Comme la déri- 

d Log.v vée (04) " 

d Log.CP 


d Log.v ? 
les dérivées 


d Log.  d Log.cit d Log,c4 


à Loc.” dd ce ne doser US 


peuvent être résolues en deux dérivées partielles, 
l’une telle que 
f'OJoUE 
ses ("6 ). 


prise à volume constant par rapport à la tempé- 
rature, l’autre telle que 4 


ë etes 

Log.v /o 
prise à température constante par rapport à la dila- 
tation (?). 4 


(?) Dans le calcul de la dilatation je prends en compte les 
Fic."8, dérivées totales des pulsations ows- et des facteurs gs-, 
c’est-à-dire : 


Ê : . d Log.ose _ 1 ess) ( Log. os) 
restent approximativement constantes (à quelques ( en de US Lors Los 


centièmes près) (1). dLog.v  Bos- 


d Log.gs- 1 [dgsr ù Log.gs- 
(?} Et à des écarts encore plus grands, on a : d Log.v  fgs: | 00 /» ù Log.v /6 
d Log. =; d_Log.T| 1 
d Log.v |  2]d Log.v 6° Manuscrit reçu le 11 avril 1961. 
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INTENSITÉS RELATIVES DES RAIES X DU SPECTRE L, 
EXCITÉ PAR BOMBARDEMENT ÉLECTRONIQUE, DES ÉLÉMENTS LOURDS. 


Par Marcez GOLDBERG, 


Institut du Radium. Laboratoire Curie, Paris. 


Résumé. — On a mesuré les intensités relatives des raies X du spectre L des éléments suivants: 
Ta, W, Re, Pt, Au, Hg, Tl, Pb, Bi, Th, U. On a travaillé avec un tube à rayons X sous une tension 
de 50 kV. Les raies sont séparées par un spectrographe à cristal courbé (Cauchois) et détectées 
par un compteur à scintillation : Nal({Tl) + photo-multiplicateur. Les résultats sont corrigés 
pour tenir compte du pouvoir réflecteur du cristal et de l’absorption dans l’anticathode. 


Abstract. — Relative intensities of L X-rays were measured for the following elements : 
Ta, W, Re, Pt, Au Hg, TI, Pb, Bi, Th, U. A 50 kV X ray tube was used. To analyse the rays a 
bent-crystal spectrometer (Cauchois) and a Nal(Tl) scintillation counter + photomultiplier were 
employed. The results are corrected for crystal-reflecting power and target absorption. 


I. Introduction. — Nous avons repris l’étude des 
intensités relatives des raies X du spectre L dans 
la région des éléments lourds pour l’intérêt qu’elle 
peut présenter dans l’interprétation de certains 
phénomènes d’interaction noyau-cortège. Dans la 
conversion interne L par exemple la connaissance 
des rapports de conversion Zx : Lar : Lux permet 
d’assigner une nature et une multipolarité à la tran- 
sition nucléaire. Expérimentalement on mesure 
Lx : Lar : Lux soit par spectroscopie 8 magnétique 
soit par spectroscopie des rayons X émis à la 
suite de l’ionisation des trois niveaux Z. Cette der- 
nière méthode est en particulier employée lorsque 
la transition nucléaire a une énergie peu supérieure 
aux énergies de liaison des électrons Z, les électrons 
de conversion sont alors peu énergiques et leur 
étude est difficile. 

Les”rayons X des sources radioactives sont en 
général étudiés par absorption sélective ou à l’aide 
d’un compteur proportionnel ;on n’arrive alors qu’à 
séparer les trois groupes de raies du spectre L : «, 
B, y. Dans quelques cas particuliers on a obtenu les 
spectres X par diffraction cristalline : conversion 
du rayonnement de 47 keV du RaD [1], spectre L 
de Pu et Np émis lors des désintégrations : 


Cm? —> ,Pu2s8 (2] 
et 
osAmM?4l —> ,4Np°57 [3] 


par exemple. Pour cela il faut une source intense, 
de plus le résultat n’est aisément interprétable que 
s’il n’y a qu’un seul rayonnement y converti, c’est 
le cas du RaD. Dans les autres cas on est conduit à 
faire des coïncidences pour isoler les rayons X 
appartenant à la conversion d’un y donné. On doit 
alors employer un compteur proportionnel. Il s’agit 
donc de pouvoir remonter de la connaissance des 


intensités relatives des groupes « : 6 : y aux inten- 
sités relatives X1 : X1r : Xrx des raies X des trois 
niveaux L. Xrest réparti dans $ et y, X1r également 
(plus n qui se trouve entre « et 8), Xrrr dans «et 8 
(plus quelques raies peu intenses /, s, t). Ces répar- 
titions sont indépendantes du mode d’excitation du 
cortège. Elles peuvent donc être étudiées par le 
moyen le plus commode, le bombardement élec- 
tronique et la spectrographie cristalline des raies. 
Un autre problème est de passer de X1 : X1r : Xux 
à Lx : Lu : Lux (ionisations relatives primaires). 
n’est pas aisé de le résoudre dans la région des Z 
qui nous intéresse car, outre les rendements de 
fluorescence de chaque niveau, il faut connaître les 
probabilités de transitions Coster-Kronig [4]. Les 
transitions Coster-Kronig consistent pour le 
spectre Z dans le passage de l’ionisation de la 
couche Zx à la couche Zxr où Lu ou de Zxr à Larx, 
l’énergie libérée servant à expulser un électron 
d’une autre couche 47.N ... Ces transitions ne sont 
énergétiquement possibles que pour certains élé- 
ments. Li — Lairr Miv par exemple ne se produit 
que Si : 
Erxz — Errn — EMjy 0 


(Euy : énergie du niveau WMiy lorsque Lin est déjà 
ionisé). 
C’est le cas pour les atomes de numéro atomique : 
ZI HO 722,75, 


Les rendements de fluorescence (w;) et les coef- 
ficients de passage Coster-Kronig (f.) sont des 
caractéristiques du réarrangement du cortège élec- 
tronique, indépendantes du mode d’excitation de 
ce dernier. La valeur des rapports Ar: Arr: Xwx 
dans un mode d’excitation donné où Lx : Lar : Lair 
(rapports des probabilités d’ionisation) sont connus 
contribue à la détermination des oo; et fi;. 
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Fic. 1. — Raïes L£ de l’uranium, 


Les mesures d’intensités relatives des raies X du 
spectre L sont peu nombreuses. Les mesures an- 
ciennes [5] ont été faites par spectroscopie à cristal 
et détection photographique ou à l’aide d’une 
chambre d’ionisation. La méthode photographique 
ne permet de comparer correctement que des raies 
de longueurs d’onde et d’intensités voisines. Vu la 
durée des mesures d’intensités, les chambres d’ioni- 
sations sont trop peu stables. Récemment Victor [6] 
a mesuré les intensités relatives des raies du 
spectre L pour les corps suivants : W, Pt, Au, Pb, 
Bi, Th et UÜ, en se servant d’un compteur propor- 
tionnel comme détecteur. 


II. Montage expérimental. — a) Les rayons X 
sont produits dans un tube démontable avec fenêtre 
de béryllium. La tension utilisée est de 50 kV. On a 
toujours des électrons dont l’énergie est plus du 
double de l’énergie d’ionisation des niveaux Z des 
éléments étudiés (Ez, (o2U) — 21,756 keV). 

On peut obtenir une intensité d’électrons allant 
jusqu’à 25 mA. Les cibles sont constituées par des 
dépôts électrolytiques au tampon (procédé Dalic) 
sauf pour Ta, W, Th et U où des pastilles ont été 
soudées après avoir été recouvertes de cuivre dans 
un four à vide. Dans tous les cas les cibles sont 
épaisses, c’est-à-dire qu’elles absorbent complè- 
tement les électrons de 50 keV. 

b) Les raies sont séparées par un spectromètre à 
cristal courbé à transmission (type Cauchois) [7], 
construit par M. Frilley [8]. C’est un cristal de 
mica de 12,5/100 mm d’épaisseur. Nous nous 
sommes servis des plans 100 dont la distance réti- 
culaire est de 2,55 À. Le rayon de courbure est de 
40 em ce qui donne une dispersion d’environ 
12,5 UX/mm. 

c) Le rayonnement est détecté par un compteur 


à scintillations INa(Tl) + photomultiplicateur 
51 AVP, avec montage en impulsions. Les impul- 
sions sont amplifiées puis passent par un sélecteur à 
un canal qui élimine en grande partie le bruit de 
fond du P. M. et les autres parasites. Elles sont 
ensuite transmises à un intégrateur relié à un enre- 
gistreur. 

Le cristal scmtillateur INa(T1) a 2 mm d’épais- 
seur, il absorbe la quasi totalité du rayonnement, 
qui a dans nos expériences une énergie toujours 
inférieure à 21 keV. Il n’y a donc pas lieu de faire 
des corrections de rendement du compteur. Avec la 
stabilité c’est un des avantages de l’emploi de ce 
mode de détection. 

Une fente portée par un chariot précède le détec- 
teur. Les diverses positions de la fente corres- 
pondent aux différentes longueurs d’onde. L’en- 
semble de l’appareillage a été conçu pour obtenir un 
enregistrement continu et automatique des spectres. 

La figure 1 donne un exemple de courbe expéri- 
mentale. 


ITT. Corrections. — Pour comparer les intensités 
des diverses raies on a à tenir compte de plusieurs 
facteurs de correction. 

a) L’absorption entre l’anticathode et le détec- 
teur. Le rayonnement traverse la fenêtre de béryl- 
lium du tube à rayons X (14/100 mm), le cristal de 
mica (12,5/100 mm en négligeant la variation de 
parcours avec l’angle de Bragg), 50 cm d’air et la 
fenêtre du compteur en béryllium (1/10 mm). 

b) Le pouvoir réflecteur du cristal est fonction 
de la longueur d’onde. Lind, West et Du Mond [9] 
ont vérifié qu’un cristal courbé se comportait 
comme un cristal mosaïque. Le pouvoir réflecteur 


esé alors proportionnel au carré de la longueur 
d’onde. 
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c) Les électrons pénètrent dans la masse du 
métal et le rayonnement est produit jusqu’à une 
certaine profondeur à partir de la surface de l’anti- 
_ cathode. Le problème des corrections d’absorption 
dans l’anticathode est étudié dans le livre de 
Compton et Allison, p. 86 [6]. Le nombre de pho- 
tons W,; à la surface de l’anticathode est donné par 
la formule : 


0 
M = AT à Elulæ)] e—we cts 0 de (1) 
% 


A est un facteur qui s’élimine lorsque l’on fait les 


TABLEAU I 
TNT NT RO TS SP 7O AUS OT TNTimES2 Pb 83 Bi 90 Th [92 U 
Lrrr My 
CA 100 100 100 100 100 100 100 100 100 100 100 
Lrx Miy 
o 12 117 AA TEA AIS 14: Lt 10,5 41,92 1190 10106 
Lrr Mr 
n 1,1 1,1 0,8 0,9 0,8 4,1 0,9 1,1 0,9 
Li Miv 
8, DS LE NL ER LL NN D -LG4 COST 6 40 CSI 
Lin Niv,v 
VAE 20:22 292926 944 71339 959. 999 -96.8 224 24 28,4 26,8 
Lx Mixx 
8, 8,6 6,7 6,1 2, 7e) 5,7 4,5 4,8 4,3 CA ss 
LT Mr 
Ba 6,5 5,3 &,h 3,7 3,6 de 4,3 k,3 LA . *#5,3 3,6 
Lair Orv,v 
Bs 0,28 0,25 0,93 151 21 259 3 2 A RAD El 
Lux Nr 
Bs 42 122 1,8 155 473 129 17 1190) 1,8 49 à Vs 
Lx Or 
8, 0,33 0,21 0,82 0,36 0,35  0,5- “0,9 - “1,1 0,42 
Li Ny1,vir 
gs 0,07 0,16 0,22 
Lx MY Ë 
Bo Di 0740 LOS 0,32 0,4 0,5 0,4 
Lx Miv 
Bo D26 4005 0,27 0,3 0,3 0,25 
Li Niv ; . 
Y1 9,6 12 12,6 14,7 13,9 13,9 137 12,8 11,6 11207 10,8 
Lx Nix ; 
Ye 1,9 1,9 1,3 1,6 1,8 0,95 1,6 1,3 1,3 1,4 1,05 
Li Nut j 2 
y: 3 2,8 2 2,2 2,3 1,5 1,8 1,6 1,7 4,95 - 1,15 
Dr Or, Out 
Y + Ya TÉMSE RE 07 1 1 1 0,9 0,7,20,6 0,77 
Lux Nr à à ; 
7 0,82: We 0:39 010:5 0,45 0,65 0,5 0,65 0,6 0,45 0,34 
18 Orv 
F Ye 00027 DS 1 1,6 0,8 1,8 1,9 1,3 2,2 9,2 
Li Or 
Ye 0,16 0,06 0,28 
* 6, + Bz. 


He Interpolé. 
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rapports d'intensité de raies d’un même niveau ; 
T probabilité de transition ; F{u) fonction d’ioni- 
sation ; w rapport de l’énergie des électrons à 
l’énergie d’ionisation du niveau. 

On prend u fonction de la pénétration des élec- 
trons : u(x). ; 

Nous avons établi cette fonction à partir des cal- 
culs de Tsien San Tsiang [10] des relations énergies- 
parcours des électrons. 

x, est tel que u(xs) = 1. 

0 angle que font les rayons X avec la surface de 
l’anticathode. 
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Nous avons pris pour fonction d’ionisation la 
fonction empirique de Webster, Hansen et Duve- 
neck [11] : 


il 
Fu) = constante X 
«|( 


nl à 


Soient deux raies 1 et 2 d’un même niveau, C; et 
C, la valeur des intégrales de (1) qui leur corres- 
pondent. Le rapport d'intensité expérimental 
N,IN, doit être corrigé par C,/C2. Dans le cas de 
raies provenant de niveaux différents on procède 
de la façon suivante : 


TABLEAU II 


DIT MP ANITER Re 78 ePt 
Mr: 
Ba 76 79 72 79 
Min 
Bs 100 100 100 400 
Miv 
se 3 4,5 5,8 
My 
Li 8 5,5 5,1 8,7 6,8 
Nu 
Yo 22 29 21 34 
Nr 
y: 35 492,5 33 47 
Orx,Oxrx 
Ya + Ya 9 12 PRO 
Mi 
* 2,4 2,4 1,8 2 
Miv 
B: 100 100 100 100 
Nr 
y: 0,7 0,84 1,1 1 
Niv 
Lu y: 21 26 28,5 32 
Or 
Ye 0,35 0,13 
Orv 
Ye 0,7 0,6 4,3 252 
Miv 
% 12 11,7 11,7 11,8 
My 
CA 100 100 100 100 
Mi 
Be 1,2 1,2 1,8 195 
Niv,v 
Lurr Bis FBe 20,2 27 21,4 29 
Nvivi 
7 0,07 
Or 
Bz 0,33 . 0,21 0,32 
Ojv,v 
Bs 0,28 0,25 0,53 4 
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Pour une raie donnée le rapport de l’intensité du 
rayonnement détecté à l’intensité du rayonnement 
produit dans la direction définie par 0 est égal à : 


: | 
. Flu) vscote dx 
0 


[rt de 


et on applique comme correction pour deux raies 1 
et 2 ne provenant pas du même niveau : 


CilC: X K:1Ks - “4 
Les intégrales ont été calculées graphiquement. | 
Nous avons travaillé avec un angle 8 — 180. 


C 
me (3) 


79 Au 80 Hg S81TI 82Pb 83Bi 90 Th QU 
37 67 95 90 95 97 
100 100 100 . 100 100 100 100 
6% 6,7 6,8 
8,5 11 11 24 
38 17 35 27 30 a1 29 
49 26 40 33 10 51 32 
21 17 20 15 14 21 
1,9 2,4 NC RAR 2,3 
100 100 100 100 100 100 400 
1,5 1,1 1,4 NOR: 0,9 
33 29 30 31 29 30,5 28 
0,6 
3,8 1,8 39 TVR ue 5,5: 5,7 
11,5 -11 10,5 . 11,9-,14,3%- 10,7 2106 
100 100 100 100 100 100 100 
1,3 1,3 1,7 1,7 48 19. 17 
25,9 22,9 26,8 224 94 28,4 26,8 
0,16 0,22 
0,36 0,35. 0,5  %09 “A 0,42 
2,1 2,3 3 se 5,7 


No 11 


A titre d’exemple nous donnons les corrections 
pour le rapport B,/«; 


Z A 78 79 82 83 92 
CORRE ET RE 070E CET SE: 


Bien que très approximatives (hypothèse du par- 
cours rectiligne des électrons, négligence de la fluo- 
rescence due au rayonnement continu, etc...) ces 
corrections améliorent les résultats. Des expériences 
faites avec différents angles montrent une concor- 
dance convenable. Les corrections sont particu- 
hèrement importantes quand on compare 2 raies se 
trouvant de part et d’autre d’une discontinuité 
d'absorption. 


IV. Résultats expérimentaux. — Dans le ta- 
bleau I nous donnons les intensités en pour cent 
de la raie «,. Il s’agit ici de rapports de nombres de 
photons, certains auteurs [5] donnent pour les 
intensités relatives de tels rapports multipliés par 
les rapports de fréquences. L’intensité d’une raie 
comprend l’intensité de la raie de diagramme plus 
l’intensité des raies satellites lorsqu'elles existent. 
Dans le spectre L les raies satellites sont liées à la 
double ionisation qui suit un passage Coster- 
Kronig mais correspondent à la même transition 
que la raie principale. Les résultats sont donnés 
corrigés comme il a été indiqué au paragraphe 3. 


© Présent travail 
A Victor [6] 
© Théor, [12] (13 


Fig. 2. — Lx Rif. 


Les tableaux II donnent les résultats niveau par 
niveau. Les intensités sont rapportées à fa, 1, « 
our Lx, Lar, Lars respectivement. Nous obtenons 
es rapports des probabilités de réarrangement qui 
sont indépendants du mode d’ionisation. 

Sur les figures nous portons avec nos valeurs, 
celles obtenues par Victor [6] ainsi que les résultats 
des calculs théoriques concernant les transitions 
L — M effectués par Massey et Burhop [12] pour 
Au et U et Laskar [13] pour W, Pt, Pb et Bi. Los 
divergences peuvent venir de l’imprécision due à 
à l’appareillage, de la façon dont on effectue les 
corrections et enfin de la manière dont on sépare 
les raies partiellement confondues comme f, et B; 
dans Pt et Au par exemple. 

La figure 2 montre la variation du rapport 84/83. 
Ces raies forment un doublet de spin. L’état initial 
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étant le même, la règle de Burger et Dorgelo prévoit 


des intensités proportionnelles au poids statistique 
(27 + 1) du niveau Mn(pir2) pour Ba et Mir (Pa/2) 


50 


LEE ON EN EE CEE ET TON 
73 74 75 78 79 80 81 82 83 30 CII, S 


F1G. 8. — Lr Y2/vs et val Ys. 


78 79 80 81 82 83 90 


nes = Tire 


T6 19 80 8 67 83 0 Rx 


Fic. 5. — Lrrr. 


pour 63. On devrait avoir donc un rapport constant 
égal à 1/2. On voit que les résultats expérimentaux 
diffèrent nettement de cette valeur, par contre ils 
se rapprochent des valeurs obtenues par le calcul 
relativiste des intensités radiatives [12], [13]. On a 
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un exemple identique de doublet avec Y2 et Ys 
(Li — Nu, Ni) (fig. 3). ue 

. Des autres figures notons que les rapports d’in- 
tensité des raies émises lors des transitions L — M 


et N varient peu, alors que les rapports Y,/B: et, 


Bs1x, sont de l’ordre de dix fois plus grands pour U 
que pour Ta. Ceci est lié au remplissage de la 
couche électronique externe Orv, Ov dans l’état 
solide. Cette couche est complète à Z = 79 dans 
l’atome libre. La variation de Y4/Ys 


(Lx a On, Our) 
Lx — N III 
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est moins prononcée, Or, Our est complète dans 
toute la région qui nous intéresse. 

Les rapports des intensités totales rayonnées par 
chacun des niveaux X1: X1r : Xrr, qui dépendent 
des probabilités d’ionisation, des rendements de 
fluorescence et des coefficients de passage Coster- 
Kronig, seront analysés dans une publication ulté- 
rieure. 

Ce’ travail a été effectué sous la direction de 
M. Frilley auquel je tiens à exprimer toute ma 
gratitude. Je remercie également Mme Blandin, 
du Laboratoire de l’Aimant Permanent, qui a 
effectué les dépôts sous vide. 


Manuscrit reçu le 22 avril 1961. 
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ERRATUM 


Dans l’article de À. Messe sur la “ Théorie de l’effet 


photoélectrique des métaux ?”, 


(J. Physique Rad., 1961, 


22, 308), la figure 3 doit être corrigée en retournant les 
courbes expérimentales de Dickey, dont l’échelle d’énergie 
doit être E au lieu de EE, — E. Ceci a d’ailleurs pour effet 
de donner un accord bien meilleur avec la théorie et de 
rapprocher énormément les résultats expérimentaux de 
Dickey et de Methfessel. Nous remercions le Dr Methfessel 
d’avoir attiré notre attention sur cette erreur. 


LE JOURNAL DE PHYSIQUE ET LE RADIUM 


PROPAGATION D’UNE ONDE ÉLECTROMAGNÉTIQUE PLANE 
DANS UN PLASMA A LA COUPURE ET AU DELA 


Par R. PAPOULAR, 


Service de Fusion Controlée, C. E. N., Fontenay aux roses. 


Résumé. — La solution complète des équations de Maxwell pour deux types de’ distributions 
de densité fait ressortir l’insuflisance des hypothèses de l’optique géométrique pour les plasmas 
minces : a) le coefficient de réflexion d’une couche de plasma n’est voisin de 1 que pour des fré- 
quences inférieures à la fréquence nominale de coupure ; b) ces fréquences sont d’autant plus basses 
que l’épaisseur du plasma est plus faible ; c) le déphasage de l’onde réfléchie dépend de l’épaisseur 
du plasma aussi bien que de la fréquence de travail. Par contre dans certains cas (au moins), le 
déphasage de l’onde transmise est pratiquement le même que celui calculé en optique géométrique. 
Une relation générale entre l’amplitude et le déphasage du champ électromagnétique est démontrée 
en annexe. 


Abstract. — The complete solution of Maxwell’s equations for two typical dersity distributions 
makes clear the shortcomings of the geometrical-optics hypotheses for thin plasmas : 

a) the reflection coefficient of a layer of plasma is very different from 1 unless the working fre- 
quency is lower than the nominal cut-off frequency and ; 

b) the more so the thinner the layer ; ÿ 

c) the phase lag of the reflected wave depends on the plasma thickness as well as on the working 
frequency. Nevertheless, in certain cases, at least, the phase lag of the transmitted wave is pra- 
tically the same as that which is deduced from geometrical optics. A general relation between 
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amplitude and phase angle of the electromagnetic field is also proved an appendix. 


I. Introduction. — L'une des méthodes les plus 
utilisées pour la mesure de la densité des gaz ionisés 
est celle de l’interférométrie hertzienne. Elle consiste 
a mesurer le déphasage subi par une onde plane 
traversant le volume ionisé ou se réfléchissant sur 
Jui. 

Pour interpréter le résultats de ces mesures, on 
considère habituellement que l’approximation de 
l’optique géométrique est valable, c’est-à-dire que 
la variation relative de l’indice de réfraction du 
milieu ionisé est faible pour un trajet égal à la 
longueur d’onde. 

Or, il n’est pas rare de rencontrer dans les labo- 
ratoires, des plasmas de quelques centimètres 
d'épaisseur. D’autre part, les densités électroniques 
courantes étant comprises entre 1012 et 1014 cm5, 
les longueurs d’ondes les plus usitées en interfé- 
rométrie sont 8 et 4 mm. Dans ces conditions, on ne 
peut admettre que la densité et l’indice soient très 
lentement variables, et il paraît utile d’analyser 
plus rigoureusement la propagation électroma- 
gnétique dans ces milieux. 

Plusieurs auteurs se sont attaqués à ce problème 
dès le début des études sur l’ionosphère. C’est ainsi 
que Hartree (1931) a donné les solutions des équa- 
tions de Maxwell et le déphasage de l’onde réfléchie 
dans certeins cas particuliers (réflexion totale) et 
s’est surtout intéressé à des variations relativement 
lentes de la densité ionique. A sa suite, Wilkes 
(1940) étudia les ondes polarisées circulairement, 
avec une densité croissant suivant une loi para- 
bolique. PET TR 

Sahàa et Raï (1937), puis Deb (1940) ont étudié 


l’amplitude du coefficient de transmission de 
couches ionisées d’épaisseur finie mais grande, au 
voisinage de la fréquence de coupure (répartitions 
de densité triangulaireet parabolique). Un travail 
qui semble mal connu est celui d’Epstein (1930) 
lequel a donné une expression du coefficient de ré- 
flexion dans le cas d’une répartition en « cloche » 
sans restriction sur les deux paramètres définissant 
cette distribution. Toutefois, il ne s’est pas inté- 
ressé aux questions de phase, ni à celles de l’ampli- 
tude des champs dans la région onisée. Tout récem- 
ment, Golant (1960) a cité des résultats obtenus par 
la méthode d’Epstein et concernant l’amplitude du 
coefficient de réflexion des plasmas minces. Enfin 
Budden (1957 et 1959) a mis au point une méthode 
de calcul matriciel numérique permettant de ré- 
soudre le problème le plus général à l’aide d’une 
machine à calculer. 

Dans le présent travail, nous nous proposons 
d'étudier les coefficients de réflexion et transmission 
en amplitude et en phase ainsi que la répartition 
des champs dans la région ionisée, et ceci pour deux 
types de distribution de densité qui se rapprochent 
assez de la réalité. 

Hypothèses : Pour faciliter les calculs, nous avons 
fait les hypothèses simplificatrices suivantes : 

1) La densité électronique et, par conséquent 
indice de réfraction ne dépendent que d’une va- 
riable d’espace x, la propagation s’eflectue le long 
de l’axe et le plasma est supposé illimité dans les 
directions y et z : le problème est donc ramené à une 
dimension. Cette hypothèse est justifiée par le fait 
que les interféromètres généralement utilisés sont 
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pratiquement des filtres d’ondes planes à cause de 
la directivité des aériens émetteurs et récepteurs. 

2) On ne tient pas compte du champ magnétique. 
— Pratiquement, il y a toujours un champ magné- 
tique perpendiculaire à x, mais le vecteur élec- 
‘ trique de l’onde hyperfréquence peut être dirigé 
parallèlement à ce champ, auquel cas ce dernier 
n'intervient pas dans la propagation. 

3) On ne tient pas compte de la dissipation d’éner- 
gie par collisions (conductibilité électrique nulle). 
Cette simplification est d’autant plus justifiée que 
le milieu est plus chaud et moins dense. 

Pour résoudre le problème posé, il est également 
nécessaire de faire une hypothèse sur la répartition 
de la densité électronique en fonction de l’abscisse x. 
Nous considérerons deux types de répartition : le 
type parabolique et le type trapézoïdal. 


II. Répartition parabolique. — Il est certain que 
cette répartition n’est pas parfaitement réalisée 
dans la pratique ; toutefois, elle ne doit pas être très 
en défaut, tout au moins dans les régions centrales 
du plasma. 

D'autre part, avec cette en la solution 
rigoureuse des équations de Maxwell ne s’exprime 
au moyen de fonctions connues et tabulées que 
dans des cas particuliers. Parmi ceux-ci, nous ne 
considérerons que celui où la densité critique ne 
correspondant à la fréquence de travail est égale à 
la densité maximum du plasma ; c’est ce que nous 
appellerons condition de coupure nominale. Nous 
prendrons, pour origine des x, l’abscisse de ce maxi- 


; 


L 
LD LE 


x, Oz, 


& L 
Fic. 1. — Dzasité et permittivité en fonction de la distance 


mum. On aura donc (voir aussi fig. 1): 
densité relative 


nfne = 1 — (xfx,)? pour — x, < x < x; 
= 0 Pourriel 


permittivité relative 
e = 1— nfne — a? [ri 
indice du plasma 
i=Ve= xx, 


indice du milieu environnant (vide ) — 1. 
L'hypothèse n(0) — ne est à peu près réalisée, 


N° 


chaque fois qu’on cherche à mesurer la fréquence 
de coupure d’un plasma donné. 


II-1. ÉQUATIONS DU PROBLÈME ET FORME DE LA 
SOLUTION. — Les champs électrique et magnétique 
sont de la forme 

8 — Ex) e0t 


Les équations de Maxwell en rotationnel devien- 
nent 


et J — H(x) efot, 


dE | 
eee ge 1 
dx 1640 H ( ) 
dH ; 
er du cE (2) 


où pu, : perméabilité absolue du vide, 
.&  permittivité absolue du vide, 
© : 2Tf, 


: fréquence de travail (= fe, ici) 
et d'où l’on tire 
d? E 
Le Pass _ SR —10 (3) 
d’H 1-4 4H. “0? 
dede (a) 


où € — vitesse de la lumière dans le vide. 


Be 
V0 Lt PSE L 2 

Ces deux équations appartiennent à la même 
famille d'équations de Bessel. 


J'+ Tu + bem y = 0. (5) 


Cette équation se ramène en effet à l’équation (3) 
si l’on fait 


y=E a = 0 
= [Cri = (ka) = (2xfae); m2 
et à l'équation (4) si l’on y fait 
y = H a = — 2 
= &?[C?af = (k]æ,)? = (2r/ax,)? m = 2. 


Or, la solution générale de (5) est de la forme 


1—a Ca RD 
= ae Bin | OU 
m+2 Lm +2 


où B est une combinaison linéaire de deux fonc- 
tions de Bessel linéairement indépendantes. Or 
dans le cas présent, l’ordre de ces fonctions est 
v = (1 — a)/(m + 2) = 1/4 ou 3/4 selon qu'il 
s'agisse de la solution de (3) ou de (4). Comme ces 
nombres ne sont pas des entiers, on peut prendre 


B, = MI, + NJ-, 


car alors les fonctions J, et J_, sont linéairement 
indépendantes, 


No 11 


Posons done Z — (/b/2) x? — argument des 


fonctions de Bessel 


E = 2% MJn(Z) — NJ (2) 
= 244 b—U8 ZUAÏMJ,,,(Z) — NJ_;n(2)] (6) 
He 2?[PJ ya (Z) Er) QI sn(2)] 


= 998 M8 ZHAIPJsn(Z) — QJ_on(Z)]. (7) 


Portons maintenant (6) et (7) dans (1), en tenant 
compte des relations 


d 
az CE Jia) = ZE sta 


d 
dZ Fe Ja) = = 


3/4 


ainsi que de 


der be ZEN AZ: 


Il vient 
PE VEN rer M 
(ON ob 
D'où 
Ta Pet À = 
Le SIN aa(Z) + MI s4(2)] 
Oo 
Ve e 
1= DIE ZEN Jsa(Z) + MJ_sn(2)1 (8) 
PT 


Tout ceci n’est valable que pour x > 0. 

Pour x < 0, posons x —— X avec X > 0. Les 
signes « moins » disparaissent dans (1) et (2) mais 
(3) et (4) ne changent pas de forme. 

Si l’on pose Z —(/b/2) X?, on trouvera pour Æ 
et H des solutions de la forme 


E = 2H VUE ZM Jp — No Ja) 


on pale ZEN PT EME Et 


qui sont de la même forme que (6) et (8), au signe 
près de A. 

Or, ces solutions doivent se raccorder aux précé- 
dentes pour 


= X = Z— 0: c'est-à-dire : 
E_(0) — E(0) (9) 
et 
H=\0) = H(0) 


Mais on sait (Jahnke-Emde-Loch, 1960, p. 134) 
que | 

(Z13)Y (ZI2)+? 
OI Tv +1). 11[{v + 2) 


Jy(Z) — 


[ 
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Il en résulte que 
Lim J;y4(2Z) = Lim J;4(2) = 0 
4 1/4(2) ds sa (2) 


91/4 
Li 1/4 PER NT 
Lim Z%4J he 
a —3/a (2) = 0! T1/4) == 0,422. 
Les densités (9) donnent alors 
MEME EN =EON: 


Finalement, on peut prendre 
E_{X) = — 204 b—U8 ZUAIMJ ,(Z) + NJ_n(Z) 
Jo18 


(10) 


MI; (Z)] 


— — 9941 (11) 


ZHAINJaa(Z) — 
CON 

Les deux constantes M et N seront déterminées 
par les conditions aux limites, c’est-à-dire, en 
l’occurence, l’onde incidente, Æ;, au plan d’abs- 
cisse — x, et la présence du vide de part et d’autre 
du plasma. 


IT-2. COEFFICIENT DE RÉFLEXION. — Pour con- 
naître le coefficient de réflexion du plasma il est 
utile de calculer d’abord «l’impédance » 

— 2) = EC 2) 2%) 
présentée par le plasma dans le plan — x:. 
Dans ce but, écrivons que l’impédance présentée 


au plasma dans le plan + x, est celle du vide, 
c’est-à-dire 


2 No le Zi Er 12 = ra JA 7; 
D'où en utilisant (6) et (8) et posant N/M =S, 
[J1a(Za) — SJ 17a(Z1)) 


Pour = 7; 


(3) ge e S J3Ja(Z1) cn J_3a(Z:) 
On entire 
SG — J_sa(Z1) + 1J 114(Z1) È 
— Japa(Zs) + JJ ra (Zi) 


En portant cette valeur dans (10) et (11) on 
obtient 


# Mecs J_;4 a J(J'i/a J'aja TT J— 1/4 J__3,4) 
2J ya Ja + J(d'aya J'aja — d—aya J_3/4) 


(12) 


z(— æ;) 


où l’argument des fonctions J est toujours Z:. 
Le coefficient de réflexion (en amplitude) s’ob- 
tient par la formule bien connue 
Ne Re Rene 
On trouve 
Jia Jia — J'apa J—3)a . 
_ J'i/a 7. + J'aya J=sja + (a Jaja — J__1ja J— 313) 


Les expressions asymptotiques de r et @ s’ob- 
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tiennent aisément à partir de celles des fonctions 
de Bessel. : 


2 " 
Pa CRAN Vos (z—° 9 


Jsn(2) 


TrZ 2 
> 17 
Je A Z +53) 


On trouve alors 
Jim ep = (1/V2) e—i2a—n/2) (14) 
1 >CO 


et 


Lim p — 0.eir/?, 
Z1—00 


F6. 2, — Coefficient de réflexion p — r e—? en amplitude. 


La figure 2 représente les variations du module r 
et de l’argument + du coefficient de réflexion. Ces 
résultats montrent que : 

1. Même lorsque la densité maximum du plasma 
atteint la valeur critique correspondant à la fré- 
quence de travail, l’onde incidente n’est pas entiè- 
rement réfléchie. La moitié de la « puissance inci- 
dente » traverse encore le plasma, même si ce dernier 
est très épais. Naturellement le coefficient de ré- 
flexion tend vers zéro en même temps que l’épais- 
seur du plasma. 

2. À un déphasage près de x]2 en avant, la 
phase @ de l'onde réfléchie par rapport à l’onde inci- 
dente a la même valeur que si la réflexion était loca- 
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lisée dans le plan médian x — 0 et si les lois de 
l'optique géométrique étaient encore valables. 

En effet dans l’approximation de l’optique géo- 
métrique, le retard de phase de l’onde réfléchie par 
rapport à l’onde incidente, dans le plan — x, est, 


(rx) (2 [ia Sr NA 
At 


Par 


-3 


Fi. 3.— Coefficient de transmission en amplitude 
TES te—1% 


La figure 2 montre que, aux erreurs de calcul 
près, la relation 
o = 22, — rf2 


est valable non seulement pour les grandes valeurs 
de x, JA mais pour toutes ses valeurs. 


II-3. COEFFICIENT DE TRANSMISSION. — Si l’on 
désigne par Æ;, l’onde progressive incidente en 
æ = — +1, par E(— Z;) celle qui résulte de la com- 
position des ondes incidente et réfléchie, en 
æ = %,, et enfin par Æ(Z,), l’onde transmise au 
delà du plasma, en &æ — x, le coefficient de trans- 
mission en amplitude est 


rte = Elx)[Ei = (E(x) [E(— 31). (E(-— «3) JE). 


Mais 
E(-— x;)[Ei = 1 + p(— x) 


par définition. Et, d’après (6) et (10) 


E(x;) J'11a(Z1) — SJ 178121) 
E(— %:) J'a(Z1) + SJ ya(Zs) 


En substituant la valeur trouvée précédemment 
pour #, on obtient finalement 


ER Lee e—il(r/2)+ Arc tg.(U/T)] 


FT vT vi (rs 


où 
P = Jin dan + Jeuya Jap 
T' = Jia Jyyn + J'oya Ja 
= Jia sa + Jin ds; 


argument des fonctions J : Z;, 
Il est facile de vérifier le principe de conservation 
de l’énergie 
TER = A, 


D'autre part, 


. re . 
et | Lim +— (1/V2)e—124 | 
Z1>00 


La figure 3 représente les variations du module t 
et de l’argument Ÿ du coefficient de transmission. 
Là encore, on constate que l'expression 


= 27: 


(16) 


du déphasage de l’onde transmise obtenue en optique 
géométrique est pratiquement valable pour toutes les 
valeurs du paramètre Z, — rx;fx, sauf peut-être 
pour z\/À < 0,5. 


II-4. VARIATIONS DU VECTEUR ÉLECTRIQUE A 
TRAVERS LE PLASMA. — Prenons comme référence 
le vecteur Æ(0), dans le plan æ — 0 


E(0) = — up 2 ms 
(3/4) 


où S'a la valeur calculée plus haut ; d’où 
E(a) _ ZW Jia — SJ 114) 


E(0) op 
l'(3/8) 
10327 (un 4) pour æ > O0 (17) 
et 


Ex) [E(0) = 1,03 ZY4 (J_1n + JinlS) pour æ < 0. (18) 
La figure 4 représente les variations de la phase 
du vecteur électrique à travers un plasma « à la 
coupure nominale », d'épaisseur 4X Pour 
æ <—2Xet x > 2X, la phase varie suivant des 
lois bien connues correspondant respectivement 
aux cas d’une onde réfléchie superposée à une onde 
plane et au cas d’une onde progressive seule. I est 
difficile de parler d’ondes transmises et réfléchies 
dans un milieu d'indice variant rapidement. Tou- 
tefois, il est remarquable que pour x > 0, la phase 
varie sensiblement comme pour une onde pro- 
gressive traversant un milieu d'indice lentement 
variable (courbe en pointillé). 

Au contraire, pour x < 0, on retrouve les oseil- 
lations de phase caractéristiques des ondes station- 
naires comme si, dans cette région, se superposaient 
une onde insidente et une onde réfléchie en æ = 0 
(ou autour de ce plan). On notera que la période 
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spatiale des « oscillations » varie, en gros, propor- 


tionnellement à l’inverse de l’indice moyen du 
plasma. ( 


Fic. 4. — Phase du vecteur électrique complexe Æ(x) dans 
un plasma parabolique, d’épaisseur 4x « à la coupure ». 
L'origine des phases est définie par le vecteur Æ(0). La 
courbe pointillé représente les variations de phase de 
l’onde incidente dans l’hypothèse de l’optique géomé- 
trique. 


CIE) 


A /'X 
— Allure des variations de lamplitude |E{(x)|, 
déduite de la relation (19). 


rc 


Pour finir disons quelques mots des variations de 
l'amplitude de E(x). On peut évidemment les dé- 
duire des relations (17) et (18). Mais on peut égale- 
ment en connaître rapidement l’allure en se basant 
sur une relation générale, démontrée en annexe, 
pour la propagation des ondes planes dans les mi- 


lieux d’indices variables. 
[E(x)|? = Cie (defdz) 1. (19) 
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Cette formule a été utilisée pour dessiner la 
figure 5, qui rappelle les variations de la densité de 
probabilité d’une particule à travers un mur de 
potentiel. 


IIT. Répartition trapézoïdale.. — Comme le 
montre la figure 1, la répartition trapézoïdale de 
bases 2x, et x,/2 est une assez bonne approximation 
de la répartition parabolique et, partant, de la dis- 
tribution la plus courante de la densité électro- 
nique. Elle a aussi l’avantage de se laisser traiter de 
façon rigoureuse et, enfin, la méthode d'analyse dont 
elle est justiciable peut s'étendre à une distibution 
quelconque, laquelle peut être toujours représentée par 
un nombre aussi grand que l’on veut de segments de 
droite, où si l’on veut, de trapèzes superposés. 


III-1. ÉQUATIONS DU PROBLÈME ET FORME DE 
LA SOLUTION. — Quelque soit le segment repré- 
sentatif de la courbe de l’indice dans un intervalle 
donné, on pourra toujours, grâce à un changement 
convenable de l’origine, se ramener au cas d’un 
segment porté par une droite passant par l’origine, 
donc d’équation 

AT: 


Les équations à résoudre deviennent 
d’E 


+ Kk2 AxE = 0 (20) 


dx? 
dE 
a = jou. (21) 


Cinq cas peuvent se présenter selon les valeurs 
de æ et À, en supposant que l’énergie se propage 
toujours dans la direction des x > 0. 


III1-1-1) 4 > 0,x > 0 (fig. 6.1). 
En se reportant à la solution de l’équation (5) et 
en posant #? À — b on obtient 


E = 21? B3[(2/3) V6 x°/?] 
= YZYS[M .J'ys(Z) es N.J_;3(2)] 
où 


y = (8/2)V8 DU = (3/2)/8kUS AUS (29) 


et 
Z = (2/3) bV? xl? = 2/3 KA? œ92, 


M et N sont des coefficients qui dépendent des 
conditions aux limites. Des manipulations ana- 
logues aux précédentes conduisent à 


H = JOIN .Joja(Z) + M. 23(2)] 
où 
5 — (3/2) D Jous = (8/2)%% FU A+ Joue (23) 
Notons que 
YI8 = (3/2)7%8 (x/4)/8 20 
où 
Zo = Vols, impédance du vide 
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Les relations (22) et (23) peuvent se mettre sous 


la forme 
FA re 


{à 4 za | PET 
H ARE AD PE 


= (x) = 2 (x) 


où l’argument des fonctions J est naturellement Z 
et où @,(Z) est la matrice de transfert du segment 


@ (@) 


LE 1Ë 


(24) 


PT 


o 
| A 


(b) 5 
(a) 2 
(b) 
Fic. 6. — Segments représentatifs de la permittivité «. 


considéré. Si les extrémités du segment considéré : 
sont désignées par les indices m et n, on aura, 


(a) = Qu x Qt x (pi). 


I11-1-2) À > 0,x < 0 (fig. 6.2). 


Posons À —— x > 0. Les équations (20) et (21) 
deviennent 
2 
D = (#4) XE = 0 
dE 
dx — ]J00t4 H 
En posant 6 —— b > 0 on voit que la solution 


pour E est une fonction de Bessel modifiée 


T+a(2) Æ jEvs J+113(12) 
où cette fois, 


Z = (2/3) VB XV? — (2/3) KA? 92, 


® Not 


Rappelons que 


d1 il 
2 r& 7 (— Liys13+ ZI 3/3) 


et 


(— L_:/313 SF Zlos). 


La matrice de transfert du segment est alors 
-donnée par 


Qa(Z) = za | 


YLis ] 
1928 Loys 7) 


92 L_23 


ITIT-1-3) À <0,x > 0 (fig. 6.3). 
La transformation 8 = — b — — K£? A > 0 
conduit à la matrice de transfert 


Ylis — Yl_}s ] 


= 1/3 
ASS (az TL; — 1528 Lys So) 


où l’argument des fonctions Z est 
Pr 2 VB: ka k|A|7/? 7/2 
3 3 


et où |A] remplace À dans les expressions (22) et 
(23) de y et à. 


III-1-4) À < 0, x < 0 (fig. 6.4). 
La transformation X. = — x = 0 conduit à la 
matrice 
Pr YT —YJ— 
Q(Z) = ZV3 ( SZ Jun C4 RER (27) 


où l’argument des fonctions J est 
D] 
Z =2R|AN? 3 


et où |A] remplace À dans les expressions (22) et 
(23) 


ITI-1-5) À — 0. 
a) e —0 (fig. 6.5a). — Les équations (20) et (21) 
donnent 
E = Mxz + N 


H = jM]ot. 


La matrice de transfert est ici 


œ 1 
Qsa = Fi ie 
oo 
b) e —B = CZ 0 (fig. 6.5b). 
Si l’on pose x? —=— k? B on trouve aisément 


E = Mer + Nec 
H = 1 (Mess Nec), 
op 


0 
La matrice de transfert devient 


, exz eT ax 
= j; j; - 2€ 
Qsb = += Lea exx WE e ax At 
OU oo 
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III-2. RACCORDEMENT DE DEUX SEGMENTS DE 
MÊME PENTE. — a) Soient deux segments des 
types (1) et (2) respectivement. Quels sont les rela- 
tions entre les coefficients M et N correspondants ? 
Les équations de raccordement (9) se traduisent 
sous forme matricielle 


(or) = tort x (0) 


Mais d’après Jahnke-Emde, p. 207, 
limzJ,. (2) = Limz/]e (7) =us 0,9% 
Fa 1/3(2) = 1/3(2) Fe 


Lim Z8 J_gulZ) = Lim Z#2 1 sja(Z) = v Æ 0,59. 


20 
On en tire 
M) _{[—1 OM, 
(3) À 0 1) (x): 
b) Si les segments sont des types (3) et (4) res- 
pectivement, on voit de la même manière que 


G)= Us 5) 
Na (0 VAN 

III-3. SEGMENTS SYMÉTRIQUES PAR RAPPORT A 
L’AXE æ — 0. — L'énergie se propageant toujours 
de — co à + co, soit UL la matrice caractéristique 


qui, lie les vecteurs Æ:, H,, en amont, et Er, Hh, 
en aval d’un segment donné : 


Ex 3 4 

Ga) = 2 Gn)- 
Une relation analogue peut être écrite pour un 
segment symétrique du premier et que nous carac- 


tériserons par un signe prime. Il est démontré, dans 
l’annexe 2, que 


1 0 


(30) 


(31) 


sx EL 0 ) 
sw #( a ects "4e Den le 62} 
Si les deux segments symétriques sont adjacents 
et si 
a b 
Ci (° d 


on trouve aisément la matrice caractéristique de 
l’ensemble : 


1 0 


nee eee) a ral 
_ fad + bc 2bc 
Rs ( 2ac ad + )" 


En tenant compte de la propriété bien connue des 
quadripoles : [U] = 1. 

Remarques : 1. Les caleuls sont simplifiés par la 
règle suivante : l'inverse d’une matrice carrée de 
rang 2 s'obtient en inversant les positions des 
termes de la diagonale principale et en changeant 
seulement les signes des termes de la 2° diagonale, 
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puis en divisant tous les termes par le déterminant 
de la matrice initiale. 

2. Les tables numériques des fonctions de Bessel 
se trouvent dans Jahnke-Emde-Lôsch (1960) ou 
dans les références citées par ces auteurs ; voir en 
particulier les tables du National Bureau of Stan- 
dards, New-York, 1948 et 1949, p. 413 et 365. 


TII-4 RÉPARTITION TRAPÉZOIDALE ISOCÈLE. — 
Ce cas est illustré par la figure 7 où l’on suppose 


F1G. 7. — Répartition trapézoïdale isocèle. 


que le vide règne en x > x, et en x < x, (e — 1). 
Cette répartition est entièrement définie par les para- 
mètres x, et B. 

La matrice caractéristique WU? est donnée par 
la relation (33), ou U est la matrice caractéristique 
de la moitié gauche (par exemple), soit 


U0 = Qu(Zn) X Qr (0) X Q:(0) x Q5 (Zo) 


X Qslaxa) X Q5 (0). 


Les matrices Q3, Q, et Q; sont données par (26), 
(27) et (29), respectivement ; d’après le $ III-2 on 
peut poser 

— 1 0 
Eee 


Qs (0) x Qa(0) 
, les arguments qui fi- 


(34) 


A(1 — B) 
ST 
gurent dans (34) sont 


Comme |A] = 


Zn=kf21—B)  Za=; Le de 


(1 — B 
aura = —V— Bkx;[4. 


On voit que tous les résultats peuvent être donnés 
en fonction des paramètres x, [x et B. 

Or la région GK est parcourue par une onde 
incidente vers la droite (e-%*) et une onde réfléchie 
vers la gauche (e*), En toute rigueur, le plasma 
n’est donc jamais parfaitement opaque ; ce phéno- 
mène est évidemment analogue à l’effet « tunnel » 
de la mécanique quantique. Toutefois, dès que 


ar1]2 = x V=Dzrily 2,3 ea /2 << 0,1 
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et l’on peut admettre que 
EclHe = Jouole = iZolV D}, 


en négligeant l’influence de la moitié droite du 
plasma et l’onde réfléchie en Æ ($ IIT-1-5b). 

Il est important de remarquer que cela ne peut 
se faire au voisinage de la fréquence de coupure 
nominale (B Z 0) pour des plasmas d'épaisseur infé- 
rieure à la longueur d'onde de travail (cf. Epstein). 

Dans ce qui suit, nous admettrons, sauf indi- 
cation expresse, que le plasma est assez épais et la 


Not 


#4 


fréquence de travail très inférieure à la fréquence. 


de coupure nominale : ce sont les conditions d’un 
plasma parfaitement réfléchissant. 


III-5. PLASMA PARFAITEMENT RÉFLÉCHISSANT. 
— Ce cas est illustré par la figure 8 où l’on suppose 


Fic. 8. — Plasma parfaitement réfléchissant." 


que le vide règne en x < æp et où l’origine des 


abscisses a été prise dans le plan où la densité 
effective est égale à la densité critique corres- 
pondant à la fréquence de travail. La matrice 
caractéristique est alors donnée par 

1 


VB QiZn) xX(T 4) X @ Ze) 


s 120 b \ 
Zn sh VE 
Sie uv ; C d 
120 Ve 
où 
a = Jiy3(Zn) Leys(Za) + J_ 152) L_ 376) 
b — J'13(Z2) T_y3(Z0) Ar J_:/3(ZD) TL:3(Z a) 
© = J_3(ZD) Lays(Za) — J'ays(ZD) 1 ys(Za) 
d = J__33(ZD) L=;y5(Za) — J'aya(ZD) Liy5(Za) 
2, 7,82 _T%D 
Zp = 3 AVIAlaD —= GET , 
_2 Taie ET LA = TET 32 TD 
Za —<kV]A| xQ de Mn re, 


l’impédance présentée par le plasma au point D est 
alors 
(a + b) 


ZD = JZ0 7 — 


(d + c) 


| N°11 
_et le coefficient de réflexion, par conséquent 


Fe Re ee (Our de Je ie nie 
pp = 1 X e—/°p C+rd+juaté 
avec 
pD = 2 arc tg an 

De nombreux expérimentateurs mesurent l’angle 
de réflexion @, (retard de l’onde réfléchie par rapport 
à l’onde incidente en D) et ils en déduisent la dis- 
tance DO = æp. Il est donc intéressant de con- 
naître la relation entre op et np. 

Quand #p/À — 0 (transition rapide entre le vide 
et le plasma) 


(35) 


Zp et Za— 0 
et J,'Z) et 1/2) see 2) 
On en déduit 
Lim op = — 2 arc tg (dfa) = — 2 arc tg V— B. 


æD/À—0 
Ainsi, pour B variant de 0 à — co, la limite de on 
varie entre O et — x. Elle est égale à — x/2 pour 


B ——1(f = f:/V2). On aboutit aux mêmes ré- 
sultats en considérant la relation bien connue 


D CG—= 01e 
qui n’est valable que pour une transition brusque 
d’un milieu à un autre. 

A l’autre extrême, quand æn/À — co, Zn et 
Za = Zp(— B)#? tendent vers l'infini (à moins 
que B — 0, coupure nominale) 

PRRUIZ) — Vanz cos © où 0 << 7 — Ê re 5) # 
Z—co 2702 


eZ 


7° V2rZ 


Il en résulte que a et b tendent vers 


e£ 
7 (cos 0,73 + cos 0,3). 


Tandis que cet d tendent vers 


#2 
TZ (cos 0,73 — cos 02). 
d’où 
Lim ep = — 2 are 18 
Dr 
(cos O3 — cos Os) 
10 teen ae TP als) 
F9 8 {cos Os + cos Os) 
=27n (37) 


Or, le retard de phase d’une onde qui, partant 
de D, se réfléchit en O sans changement de 
phase puis revient en D est 

TE RO 


AE 


Une 
r JD 42 x SE VIAÏ +" — 22. 


PROPAGATION D’UNE ONDE ÉLECTROMAGNÉTIQUE PLANE 
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Ce résultat est à rapprocher de celui du $ II-2, 
équation (14). 

Le raccordement entre les deux valeurs extrêmes 
(36) et (37) de œp se fait par une courbe sinueuse 
voir la fig. 9). 


# 
450° # 
74 
2 
400 A r 


-200° 


F1G. 9. — Retard de phase de l’onde réfléchie 
sur un plasma parfaitement réfléchissant. B — emin 


Cette courbe tend vers l’asymptote (37) d’autant 
plus rapidement que |B| est plus grand : la diffé- 
rence entre la valeur réelle et la valeur asympto- 
tique devient négligeable dès que Zn et Za dépas- 
sent la valeur 2. Si — B < 1 il faut pour cela que 


1 


ap [À > = c ren TEE (38) 
Par exemple, si fe/f — 1,2, B = — 0,44 et la 
condition (38) devient 
ap [A > 1,6. 


On voit que, pour des fréquences peu inférieures 
à fc, le raccordement ne se fera que pour de grandes 
valeurs de æp/À. A la limite pour f = fe, B = 0 le 
calcul montre que pp tend asymptotiquement vers 


2ZD — 7x6 (39) 


Ainsi, pour — 1 < B < 0, ® se situe grossiè- 
rement entre les valeurs (37) et (39) et ne dépend 
pas seulement de xp}. 


ANNEXE 1 


Relation entre la phase et l’amplitude d’une onde 
plane. — La solution de l’équation (13) de propa- 
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gation du champ électrique peut se mettre sous la 
forme 
E(x) —= E, eD(x) 
où ® satisfait à la condition 
D" + D? + K? efx) — 0. 


Posons D’ = «& + 76 ; l’équation précédente se 
subdivise en 
B’ + 2a8 = 0 
et œ + a2— PR? + ke — 0 


L'intégration de la première ce ces deux équa- 
tions donne 


Log (8/80) = — 2 f «dx 


8 e+2%/adx — @, — Cte d'intégration 
Mais 
d 
Sr [äm(D)] — rs 


où © — Jm{[®] est la phase de Æ au sens habituel 
du terme. 
D’autre part, 


Î « dx — Red) 


et, par conséquent 
etats — |E|?]E3 
d’où 


lE|2.do/dx = Bo E9 (19°) 


qui équivaut à (19). On peut généraliser pour un 
milieu dissipatif en écrivant & — & — Je; on obtient 
alors 


Er = 888 + fr er dx (19) 


Notons que 
E? doJdx = — wuo RelE X H*). 


Re(E X H*) = puissance transmise par unité de 
surface à travers le plan d’abscisse +. 


JR aff de = 
(eo) 


sommée dans un cylindre de longueur x et d’aire de 
section droite àgale à 1. Il en résulte que 


Up © ÏX puissance con- 


EË = — ou, X 


[puissance transmise en x + puissance consommée 
de 0 à x]. 

La relation (19) r’exprime donc rien d'autre que le 
principe de conservation de l’énergie. 

Remarques : 1. ® est la phase du champ élec- 
trique résultant et non celle d’une quelconque onde 
incidente ou réfléchie. 
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2. La relation ci-dessus démontrée serait évi- 
demment valable pour Æ si e était constante et 
variable en fonction de x. 

3. La constante d'intégration B, est caractéris- 
tique non seulement des propriétés du milieu mais. 
aussi des conditions aux limites. 


ANNEXE 2 ‘ 


Matrices de transfert de tranches de plasma sy- 
métriques par rapport au plan central. — Consi- 
dérons la tranche de plasma comme un quadripôle 
de matrice W (fig. 10). La matrice U', d’une 
tranche symétrique par rapport au plan central 
(fig. 10) est identique à celle de la tranche initiale 


E, ,H, £, 4, £,4, E,4 
(b) (e) 
F1G. 10. — Quadripôles représentatifs de tranches de 


plasma symétriques par rapport au plan central. 


lorsque le sens de propagation de l’énergie est in- 
versée. Or pour inverser le sens de propagation, il 
suffit de changer partout H en — A (fig. 10). Or 


: (x (2) 
G)-6 104): 6-6 CA 


On en déduit 


A ee Vi Fe E, 
(; a) a (, ET 4 > au (2 D: 
La matrice cherchée est donc 


ne ne LT & mn 


qui est équivalente à celle définie par (32) 


AL’ = ? 


Manuscrit reçu le 20 mai 1961. 
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CALCUL DES INTÉGRALES DE RECOUVREMENT 
DES FONCTIONS D'ONDE DE DEUX OSCILLATEURS HARMONIQUES DIFFÉRENTS 
APPLICATION AUX INTENSITÉS DANS LES SPECTRES DES MOLÉCULES DIATOMIQUES 


I 


Par JEAN RUAMES, 


Laboratoire de Physique Expérimentale de la Faculté des Sciences de Bordeaux (*). 


+ 
Résumé. — Étude des éléments de matrice < vifrrlus > — L pr) re Palr) dr quand 
—0C0 


War) et War) sont des fonctions d’onde de deux oscillateurs harmoniques décalés et de fré- 
quences w, et &w, différentes. Relations de récurrence. Cas où n — 0 : expression simple du résul- 
tat quand ©: — &,, développement en série quand w, Æ ©,, application au calcul des intensités 
dans les spectres de molécules diatomiques (en particulier dans la séquence principale). Un deu- 
xième article traitera les cas où n £ 0. É 


+ 
Abstract. — Study of the matrix elements < v,fr#[v, > — [ pr) 70. Wir) dr when Y,, 
—C9O 


and Ÿ,,(r) are wave functions of two oscillators,out of phase, and of different frequencies w, and «2. 
Recurrence relations. Case where n — 0. Simple expression of the result, when «> = 1. 
Development in a series, when w, Æ &2 Application to the calculus of the intensities, in diatomic 
molecular spectra (specially in the main sequence). 


A second article will deal with the cases where n Æ 0. 


I. Introduction. — Dans le spectre électronique 
d’une molécule diatomique l’intensité d’une bande 
’:, 0, est proportionnelle au carré de l’élément de 
matrice 


< iles > = [7 Fat) Rat) Var) dr 


où Ÿ,, et Ÿ,, sont les fonctions d’onde de vibra- 
tion et À,(r) le moment de transition électronique 
calculé pour chaque distance internucléaire 7 en 
supposant les noyaux fixes. En première approxi- 
mation on suppose A.(r) constant, ce qui permet 
de se limiter au calcul de 


CO 
<> = if Palr) Pulr) dr. 


On peut aussi essayer de développer A(r) sous 
forme d’un polynôme en r pour mieux représenter 
les résultats expérimentaux [1]. Le calcul exige 
alors l’évaluation des éléments de matrice 


< vifrlva >. 


Dans cet article nous nous limiterons au cas 
où Ÿ,, et Y,, sont des fonctions d’onde d’oscilla- 
teur harmonique (1). Diverses méthodes [3],[4]per- 


(*) Une grande partie des résultats présentés dans ce 
travail a été obtenue par l’auteur lorsqu'il se trouvait à 
l’Institut de Physique de la Faculté des Sciences de Lille. 

(:) Dans ce cas, le remplacement classique des intégrales 
de 0 à + © par des intégrales de — co à + co se justifie 
par des considérations analogues à celles de Ter Haar [2] 
pour le potentiel de Morse. 


mettent d’ailleurs d'introduire des corrections d’an- 
harmonicité dans ces calculs. 

Le premier calcul des < v,fv, > a été fait par 
E. Hutchisson [5] en fonction de r,;, — r, et de 
x —V/w2/w,. Mais les expressions obtenues se 
compliquent vite quand v, + v, croît, surtout dans 
le voisinage de la séquence principale v, — v,. Ces 
formules ont été tabulées depuis par Bates [4] dont 
les tables à double entrée (Ag et w»/w,) sont 
extrêmement commodes. Mais elles n’existent que 
pour ?, + v, < 4 ainsi que pour 5-0, 6-0, 0-5, 0-6. 

C. Manneback [6] étudiant le même problème a 
pu, grâce à l’adoption de deux variables plus symé- 
triques par rapport aux deux états 1 et 2, recon- 
naître dans les < v.,|v, > des polynômes d’Hermite 
à deux variables et utiliser les relations de récur- 
rence qui les lient [7]. Ces relations de récurrence 
permettent de calculer de proche en proche le ta- 
bleau des < v,]v, > relatif à une molécule déter- 
minée. D’autre part G. Manneback et A. Rah- 
man [8] ont indiqué une méthode d’obtention du 
tableau des < v,friv, > à partir du tableau des 
< Dilvs >, 

P. A. Fraser, W. R. Jarmain et R. W. Nicholls [1] 
ont proposé une méthode approchée d’évaluation 
des < v,fr"[v, > à partir des < vf, >. 

Enfin l’auteur [9], [10], [11] a indiqué le prin- 
cipe d’une méthode de calcul des << v,fv, > à 
l’aide d’un développement en série. C’est cette 
méthode, particulièrement commode dans la sé- 
quence principale, que je vais exposer dans cet 
article. Des relations de récurrence généralisant 
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celles de Manneback [6] ont également été obtenues 
au cours des calculs. Un autre article sera consacré 
au calcul des quantités < v;fr"}v, > 


II. Notations. — Les deux états électroniques 
envisagés seront notés 1 et 2. Comme les deux 
états jouent le même rôle nous pourrons, comme 
Bates [4], supposer que la numérotation est telle 
que &1 < @», ©1 et w, étant les fréquences de 
vibration en kaysers des états 1 et 2. Nous appel- 
lerons r, et r, les distances internucléaires d’équi- 
libre et nous introduirons l’abscisse x =r —r, 
comptée à partir du point O0 qui correspond à la 
distance internucléaire r,. Nous poserons 


_ Vorlo;—1 
Vos/o +1 
A =(1—a°)(r, —r;) (1) 
r = Pr a 


Enfin nous appellerons « unité © » l’unité de 


longueur = V/hluwc (?) qui permet d’écrire les 


fonctions d’onde normées de l’oscillateur harmo- 
nique de fréquence w sous la forme réduite : 


ES NEA 
oule) = / ee Hate) En 0 


L’ «unité «w » que nous emploierons correspond 
à w donné par : 


Vo = (Vo; +Va.)/2. (3) 


C’est avec cette unité de longueur que r;,r2, ro 
et À doivent être évalués. Ag utilisé par Bates [4] 
est lié à À par À = Ag(i — a). 

Dans ces conditions, les fonctions d’onde Y,,(x) 
et W, (x) normées par rapport ‘à la variable x 
s’écrivent : 


= Vi — Pu(X 1) 
= V1 + € Ov (4°) 


Ta(x) 
Er (x ) 


X, = x-- {ax + A/2) 
Ya = 2 + (ar + AJ) 
(4) 


Les éléments de matrice qui nous intéressent 
deviennent alors : | 


nes = he 


avec 


2% Qu(A) Puf) dr. (5) 
Comme ®,(X,) p,(X2) ne dépend de a et de A 


6 } Cette unité s’évalue en À par 1/0,1722 Vu osi © 
est en kaysers et si 14, masse réduite de la molécule, est 
évaluée en unité de masse atomique (échelle physique). 

(8) Nous adoptons la définition habituelle des polynômes 
Hilx) = (—1)jre* ee (e—#). Il importe de 
noter que Bates [4] supprime le facteur (—1)° ce qui 


change le signe de < v,|v, > quand v,et v,sont de parités 
différentes. 


d’'Hermite : 
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No 
que par l'intermédiaire de À + 2ax. on peut écrire : 
Pn(X1) Pn(X2) = F(x, À + 2ax). (6) 


Il existe donc une relation entre les dérivées des 
éléments de matrice par rapport à a et par rapport . 
à À. Nous introduirons la notation : 


dr < vfærfu, > 


< vilzr|v, 0) = 2» dA? (7) 


et on vérifie facilement que : 


PE | Dilr le > 1 
< vilanto|v, >) o ee — a? (8) 


VA — a? 


Enfin nous appellerons <ov,jx"jv, >, et 
< v,[2"1v, 6 les fonctions de À que l’on obtient 
en posant a —0 dans les quantités que nous 
venons de définir. 


ITT. Relations générales de récurrence. — Cal- 
culons < v,|X,x"}v, > d’une part en rempla- 
çant X, par sa valeur (équation (4)), d’autre part 
en utilisant la relation (cf. par exemple [12]) : 

X qo(X [Vo Pe—1lX) + Vo +1 Do+1(X)]. (9) 

Traitons de même < v,|X, x"|[v, >. Nous obte- 
nons ainsi deux relations qui, par soustraction et 
addition membre à membre, se transforment en : 


<v,ler+ilv,> 


= ton D — Aaron, > + Vos < vilar|v, —1 > 
2V2 
HVD ET <v,+ihrln, > +Vn, F1 < vlzlv, +1 >] 


.. (10) 
2a < viantib, > + A <ul> 


1 _— 
— a LVv < vilerlvs — 1 > — Vo, < D, —1/ærlv, > 


HVos HT vlan +1> — Vas +1 <v +1llv,> 
un 


Une troisième relation s’obtient en remarquant 


que : 
+co qd 


Zétdr 


(Pol) Fo(x) 27] dx = 0. 


Il suffit d’effectuer le calcul de la dérivée et 
d'utiliser la relation [12] 


pUX) = — a [Vo Na, RE) Vo +1 Po+1lX)] (12) 


pour arr iver à: 


1+a 


[Vos < miens —1> Vo, +1 < vien, +1 3] 
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— Afevs > 
— Vo, +1< 0 + Afro, >] 
+. À Den iln, > = 0, (13) 


A œR a 
v2 


Enfin en calculant FT vlæ"[v, > à l’aide de 


la relation (12) puis en dérivant p fois les deux 
membres par rapport à À, on arrive à : 


<4 Dile|ve > +1 


ns 7 [Vos < viletl0s — 1 > 


Æ Vo, < D, — Afro, >) 
+ Vo, +1 < 0, + 1frrr, > 
— Vos +1 < vifrrln, + 1 ©), (14) 


Parmi les innombrables combinaisons de ces 
quatre relations notons en particulier les suivantes : 
éliminons entre les relations (10), (11) et (13) : 


vie luc et vins + 1: 
(1 + a?) < vifærtiu, > 


s ARS 
_ "Vo, < vlan, —1 > 
Re 2 1 


EVE < 4 — Alam > 
ñn aÂ 
SE 3 < vfer— tif > EUCE < vifztfvs > (15) 
tn rep et nv 29, 1: 
| Vo +1< 0, + 1frlu, > 
= en L Va < Dilær|vs — -1> 
si Aa 3 EN s 
- y) 10e RES [D > 
TR Li 
pe ce LS AE [De > 
D va < v, —1fælvs >. (16) 
30 < vifatilo, > et < 0, + 1]r fn > 
Vue + 1 < ville, +1 > 
ee VD, < 0, — 1m > 
1 — a À 
HS <a > 
;. n 1+ a A 
3 Dre. 4 RES 
2 a Lu) 
de ras Pi < vale —1 > (17) 
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. Enfin en dérivant p + 1 fois par rapport à A 
la relation (10) et en utilisant la relation (14) on 
arrive. à : 


< vfentilo, +0 


=? LV +1 (Us + 1) < 031 + Afarlu, — 1 ©) 


— 2 Voi(v, s 1) < 1 — Aero, + 1 > 


a VUalDs — 1) < vilztlos — 2 D 


— Vofv, —1) < 0, — 2/0, >@ 
+ Vos +1) (01 +2) < 0, + 2frrv, >) 


— Vis +10 +2)<nlm +250] (18) 

En faisant n — 0 dans les relations (16) et (17) 
on retrouve les relations TI et II de C. Manne- 
back [6]. Dans les mêmes conditions la relation (15) 
donne la relation utilisée par C. Manneback et 
A. Rahman [8] ({). Les relations que nous venons 
d’obtenir (en particulier (15)) permettent donc 
d’étendre les méthodes proposées par Manneback 
au calcul des < v,|x"|v, > pour n quelconque. 


IV. Éléments de matrice < v,[v, >.— 10 CaL- 
CUL DES ÉLÉMENTS DE MATRICE < V9 >>9 (CAS 
PARTICULIER &@y = 3). 

J’ai pu montrer [9] que, dans ce cas, l’élément de 
matrice prenait la forme simple suivante (si 
De Z V1) : 


<Dal0e > 0 = Vol fod es (A/V3) A LEA (A2/2) (19) 
Le est un polynôme de Laguerre : 


1 dr 
— —etx 4 ——{e tanta) 
x) = je" ee (e ) 


Si v, < v, on utilise la relation : 


< Dale Do = (<< Del > 0. (20) 
Je rappelle que les tables de e—#? L{(t) calculées 
par Tricomi [13] pour nr < 10 fournissent direc- 
tement < nfn >, en posant t{ — A?/2. Des tables 
des << v,[0, >, pour les séquences voisines de la 
séquence principale ont été calculées [14]. 


20 CALCUL DES ÉLÉMENTS DE MATRICE << [Ve > 
(CAS GÉNÉRAL © Go). 

Nous chercherons à profiter du résultat précédent 
et pour cela nous développerons 


< Vifve > 
V1 — a? 


(4) Rappelons la correction qu’il faut apporter à cette 
relation : Yn’et Zn” doivent être remplacés respectivement 


par Y Vn'[2et par Z Vn"f2. 


G(A, a) — 
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en série par rapport à a. En utilisant la US (3) 
on montre facilement que : 


G(A, a) = < di >0 +2 Lea 


; a? [Us 
me + cs a(21 
+. +40) #2 041 
en posant 
(ee) = < mir >P. (22) 
V1/9 > 


La relation (18) (en posant nr = p) sert alors de 


> # D Lt 
relation de récurrence pour calculer | ÿ à par- 
p+1 


tir du tableau des ( 2 . En particulier, après quel 


1/ ? 
ques transformations utilisant les propriétés des 
< vlv, >. les deux premiers coefficients peuvent 
1102 0 P 
être mis sous la forme suivante : 


Ce). = Va(0 + 1) Lu +1 —15>0 : 


UAT 


— Vo(n, + Av 103 +410 (0 (23) 


(x: — — [(20: Vs + Vs + Ve) < Dil0e 0 


Vi 


+ V{ns +1) (v2 +1) < o + 
— Voi ds <v,—1[0 —1>, 
— Vo (v, — 
— Vr,(v, — 


Alve +1 > 


1) (D1 +1) (D: + 2) LV, +2[v;—2 Do 


1) (De +1) (02 F2) <D1—2/03+2 0]. 
(24) 


La relation (18) permet aussi de trouver des 
bornes supérieures pour les divers coefficients du 
développement en série. Il suffit d’utiliser l’inégalité 
de Schwarz [<vlv, >| <° 1 et l'inégalité 


Vab < (a + b)/2. On obtient pour és) une 
p 


1 
borne supérieure de la forme f,(ÆX) (avec 


K = v, + v, + 1) satisfaisant à la formule de ré- 
currence suivante 


fotal&) = 7 PKfo(A) + (E — 2) folK — 2) 


+ (K +2) f)(K +2)]. (25) 


On trouve ainsi : 


f(K) = K 

falK) = Æ° +2 

fa(Æ) = Æ° + 8K 

fatK) = KA + 20K? + 24 
etc. 


(5) Dans la référence [10] ce coefficient était affecté d’un 
signe erroné. 
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La borne so du terme en æ s'obtient. 
ensuite par — = HE ) et le tableau 1 donne, quand 


de 04 Fi cue — 1,5), les bornes supérieures des 4 
premiers termes de la série en fonction de. 
D, +02 = K — 1. Ce tableau permet de trouver 


Jes bornes supérieures correspondant à n’importe 


quel a, par simple multiplication par (a/0,1}. 


TABLEAU 1 
BoRNES SUPÉRIEURES DES TERMES EN &@, @?, dŸ ET a‘ 
pour a = 0,1 
DIEU 1er 28 3e &e 
ORDRE ORDRE ORDRE ORDRE 
0 0,10 0,015  0,0015  0,0002 
2 0,30 0,055  0,0085  0,0012 
k 0,50 0,135 0,028 0,0048 
6 0,70 0,26 0,067 0,0142 
8 0,90 0,42 0,133 0,034 
10 4,10 0,62 0,24 0,071 
12 1,30 0,86 0,39 0,133 
14 1,50 113 0,58 0,23 
16 1,70 1,45 0,84 0,37 
18 1,90 1,81 4547 0,57 
20 2,10 2,21 1,57 0,85 
30 Discussion. — On peut faire les remarques 


suivantes sur la série (21). 

19 Cette série a pour rayon de convergence 1, ce 
qui entraîne sa validité quel que soit le rapport 
w2/w.,. En effet, en remplaçant dans les formules 
d’Hutchisson [5] « et Ax par leurs valeurs en fonc- 
tion de a et de À (x = Ÿ{(1 + a)J(1 — a) et 
Ag = Ag = A/1 — à), on s'aperçoit que la fonc- 
tion obtenue a pour seuls points singuliers relatifs 
à la variable a les pôles + 1 et — 1. 

20 Si v, — v, = v, les termes de rang p impair 
sont identiquement nuls. En effet dans ce 
cas F(x, A) (cf. équation (6)) et ses dérivées par 
rapport à À sont des fonctions paires de x. 

30 Lorsqu'on décide du nombre de termes à uti- 
üiser, 1l ne faut pas oublier que les bornes supé- 
vieures chtenues ci-dessus sont nettement pessi- 
mistes. Par exemple, dans le cas de < 2/2 À 
(cf. [10]), le terme en a? n’atteint la moitié de sa 
borne supérieure que pour À — 0et il décroît vite 
dès que A s’accroît. Des bornes supérieures plus 
strictes seraient souhaitables, même si elles devaient 
faire intervenir A. 

Ces remarques et l’examen du tableau 1 mon- 
trent que si v, — v, on a très facilement une 
bonne PPROARANOSS En ne prenant que 
< Vive 0, la borne supérieure de l'erreur (prati- 
quement égale à celle du terme en a?) n’atteint 0,05 
que pour D, = v, = 3 81 w,/o3 < 1,22. Dans les 
mêmes conditions, on peut aller jusqu’à vs = v, =6 


0,394 


Han 


| 1,581 1,732 


2,449 
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4,000 1,096 1,183 


1,874 2,000 


im 


2,646 2,828 3,000 


No 11 
0,632 0,775 
VE OMS 20710 
| 1,265 | 1,342 | 514 
2,121 2,236 2,345 
\ (IE 


Intensités dans la séquence principale en fonction de À (Sur chaque schéma les traits verticaux ont une 


longueur proportionnelle au carré des éléments de matrice < v|v >.) 


si ©2/&@3 < 1,10. Si l’on souhaite plus de précision, 
il suffit dans la plupart des cas d’évaluer le terme 
en a? (c’est alors le terme en a qui limite la pré- 
cision). Par contre, si v, = v,, il est presque tou- 
jours indispensable de calculer le terme en a. De 
plus, même si on calcule aussi le terme en a?, c’est 
le terme en a et non celui en af qui limite la pré- 
CIsion. 

La méthode que nous venons de décrire est donc 
spécialement recommandable pour l’étude de la 
séquence principale d’un système de bandes. En 
effet nous venons de voir que, dans d’assez larges 
limites A est pratiquement le seul paramètre 
réglant la répartition d’intensité, a n’intervenant 
que très peu. La répartition d'intensité dans la 
séquence principale peut donc, tant que a n’est pas 
trop grand, être représentée par une unique série 
de schémas correspondant chacun à une valeur de À 
(figure 1 établie à l’aide des tables de Tricomi). 
Cette figure permet par simple comparaison avec 


un spectre d'évaluer approximativement le A cor- 
respondant (en tenant compte évidemment du fac- 
teur de Boltzman relatif à chaque niveau de vibra- 
tion de l’état supérieur). 

D'autre part cette méthode de calcul doit per- 
mettre la construction de tables des < v,fv, > 
en fonction de À et de a d’une manière plus com- 
mode qu'avec les formules d’Hutchisson. En effet 
la complication de ces dernières les rend vite inuti- 
lisables quand v, + v, augmente, surtout dans la 
séquence principale. Ces tables seraient d’ailleurs 
plus pratiques que celles de Bates, car le choix de 
notre À réduirait beaucoup les variations des 
< v,]v,> en fonction de 4, ce qui faciliterait les 
interpolations. 

Mentionnons enfin que ces calculs ont trouvé 
d’autres applications [15]. Ils fournissent en effet 
les << v,|v, > sous une forme relativement ma- 
niable. 

Manuscrit reçu le 9 mai 1961. 
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Résumé. — Le ralentissement dans un solide d’un atome primaire animé d’une grande vitesse, 
se fait par une succession de collisions pouvant déplacer des atomes « secondaires ». On a pu 
montrer que, dans certaines conditions, il pouvait y avoir « canalisation » des chocs de long d’une 
rangée d’atomes ne donnant naissance à aucun nouvel atome « déplacé ». Ce phénomène a été 
étudié théoriquement et mis en évidence expérimentalement. 


Abstract. — In a solid, the slowing down of a high speed primary atom involves a succession 


of collisions which may displace secondary atoms. 


It has been shown that, under special condi- 


tions, a focusing process can take place along a row of atoms, without producing of a new displa- 
ced atom. This phenomenon has been theoritically investigated and experimentally demonstra- 


ted. 


Introduction. — L'effet des radiations sur les solides 
est un problème extrêmement complexe. Le compor- 
tement des solides dépend de leur nature (cristaux 
ioniques, métalliques, moléculaires ou cristaux de 
valence) et aussi de la nature de la radiation utilisée. 
Un des effets des radiations est de déplacer un certain 
nombre d’atomes (atomes primaires) qui sont lancés 
à travers le cristal avec une certaine énergie cinétique. 
Par une succession de collisions, ils vont à léur tour 
expulser d’autres atomes appelés «secondaires » hors 
de leur position d'équilibre et la structure du cristal 
va se trouver finalement perturbée par un ensemble 
de lacunes et d’interstitiels. 

Le nombre de défauts ainsi produits dans les solides 
a été évalué, ces dernières années, à l’aide de diverses 
approximations [1]. La première étape du calcul est 
l’estimation du nombre d’atomes primaires m, pro- 
duits par unité de volume au cours d’une irradiation. 
On a : 

np — Din, 6a 


où ® — nombre de particules incidentes par unité de 
| surface et par unité de temps ; 
t — durée de l’irradiation ; 
n, — nombre d’atomes par unité de volume ; 
oa — section efficace, par atome, pour les collisions 
introduisant des déplacements permanents ; 
elle dépend du type d'interaction utilisée. 


Dans une seconde étape on évalue le nombre 
d’atomes déplacés par un atome primaire. Un premier 
type d’approximation [2, 3, 4] consiste à analyser le 
phénomène comme une cascade de collisions à deux 
corps, entre atomes en mouvement et atomes au repos. 
Leur liaison avec la maille est simplement schématisée 
par l’existence d’un seuil énergétique Ea (= 25 eV) : 


(*) Travail effectué pendant un stage au C. E. N. de 
Saclay. 
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pour mettre l’atome en mouvement il suffit de lui 
communiquer une énergie supérieure à Æa ; on consi- 
dère qu’il peut ensuite se mouvoir librement jusqu’au 
prochain choc. Suivant ce modèle on peut s’attendre à 
ce que le dommage résultant soit un ensemble de 
lacunes et d’interstitiels en nombre égal (paires de 
Frenkel) répartis au hasâfd, dans un domaine relati- 
vement petit. Les mesures de résistivité à basse tempé- 
rature conduisent cependant à admettre que le nombre 
d’atomes « déplacés » (c’est-à-dire le nombre de paires 
lacune-interstitiel) est inférieur à ce que prévoit la 
théorie (= facteur 5). L’étude du désordre provoqué 
dans certains alliages par irradiation fournit une justi- 
fication supplémentaire à l'introduction des colhsions 
de remplacement dans lesquelles un plus grand nombre 
d’atomes en mouvement vient remplacer les atomes 
de la maille sans produire de défauts permanents 
supplémentaires. 

Un deuxième type de modèles fait intervenir des 
collisions à plusieurs corps.[Le petit domaine entou- 
rant l’atome frappé se comporte comme s’il était brus- 
quement chauffé et son refroidissement obéit aux lois 
de la conduction thermique dans un milieu continu 
lorsque l’excitation est suffisamment faible pour ne 
déplacer que quelques atomes/; c’est la théorie des 
« points chauds » (thermal spikes). Au contraire, 
Brinkman [6] estime que le ralentissement des atomes 
primaires de grande énergie (= 10 keV) est très rapide 
et provoque dans un petit domaine une très grande 
densité de déplacements ; c’est la théorie des « pointes 
de déplacements » (displacement spikes). 

Ces différents modèles sont difficiles à harmoniser 
entre eux et, de plus, on peut leur reprocher de ne pas 
tenir compte de la présence d’une structure périodique 
qui pourrait introduire des corrélations entre les colli- 
sions successives. Or en 1956 Wehner [7], en bom- 
bardant des monocristaux avec des ions Hg+ d’énergie 


. 
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variant de 50 à 150 eV a observé une éjection préfé- 
rentielle des atomes selon les rangées de plus forte 
densité, c’est-à-dire (110) pour les réseaux cubiques à 
faces centrées (Ag, Cu) et (111) pour les réseaux cu- 
biques centrés (W, Mo, «-Fe). Il suggère alors que la 
cascade de chocs qui prend naissance au point d'impact 
se propage préférentiellement selon les rangées les 
plus compactes. Silsbee [8] a donné le nom de « focu- 
sing collisions » à de tels phénomènes en leur donnant 
une justification théorique. Depuis, un ensemble d’ex- 
périences et d’hypothèses sur ce processus semble 
paorer une théorie plus cohérente des chocs dans les 
solides. 


Canalisation des chocs le long d’une rangée d’atomes. 
— Les premiers travaux sur la «propagation directive » 
des chocs interatomiques dans les solides sont dus à 
Silsbee. Il considère une rangée de sphères dures, de 
diamètre À, la distance entre les centres, à l’équilibre, 
étant d (distance interatomique). Si la première sphère 
est lancée avec une énergie Æ, dans une direction 
faisant un angle 0, avec la rangée, on va assister à une 
succession de collisions à deux corps, se propageant le 


FIG." 41° 


long de la rangée (fig. 1). Si 0,4, et E,+, sont les 
paramètres liés à la (7 + 1)ème sphère, on a les 
relations : 


sin:0, +: 


“CE — « COS 0,— (4 — x? sin? 0,)7/? où « — d/R (1) 


(En+ 1) [En = cos? (0% ia On+1)- 


— Le modèle choisi utilise deux principales approxi- 


mations : 

a) Les collisions à deux corps : Dans les solides où les 
distances interatomiques correspondent à un équilibre 
entre des forces de cohésion variant lentement avec la 
distance r entre les atomes et des forces de répulsion 
(dues au recouvrement des couches électroniques in- 
ternes) variant très rapidement avec r, on peut estimer 
que cette hypothèse est correcte. C’est le cas, par 
exemple, des métaux et des cristaux ioniques. 

b) Les sphères dures : Avec la forme de potentiel 
utilisée (V(r) — A er) on démontre que l’approxi- 
mation des sphères dures est valable lorsque le trans- 
fert d'énergie pendant le choc se fait avec une grande 
efficacité (en réalité il faut que T JE, > 2/8R où T 
est l’énergie transmise et Æ, l'énergie incidente). Le 
diamètre de la sphère dure est alors donné par 
V(R) = 0,5 T. Par contre si TE, est grand, on ne 
peut plus supposer que l’atome frappé reste immobile 
pendant la durée du choc. On peut estimer grossiè- 
rement que ATJT = (T]E, Bb)? où b est le para- 
mètre d'impact. 
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Ces deux approximations se contrebalancent donc 
mais non dans le cas d’une propagation linéaire des 
chocs où l'efficacité du transfert d'énergie est grande 
mais où b — 0. 

Reprenons maintenant les équations (1). On voit 
que si « < 2 cos 0,, 0,41 < 0, et l'angle 0 va décroître 
au cours des collisions successives donnant ainsi nais- 
sance au phénomène de « canalisation des chocs ». Au 
bout d’un certain nombre de collisions, nous aurons le 
transfert d’une énergie constante. 

La condition d/R < 2 cos 0, montre que l’effet aura 
le plus de chance de se produire sur les rangées denses. 
Dans un cristal cubique à faces centrées, le plus grand 
angle 0 de la vitesse initiale avec une rangée (110) est 
de 240 environ. Donc si « < 1,6 une fraction appré- 
ciable de l’énergie initiale sera canalisée le long de la 
rangée (110). 

Utilisant un potentiel de la forme 


V{r) = 0,038 exp [— 17,2 (r — d)]d] 


1 
et V(R) = > Æ,, Silsbee [9] a trouvé que dans le cas du 
2 ] q 


cuivre, x = 2 si E, = 200 eV, alors que « = 1,6 corres- 
pond à Æ£, = 50 eV (R croît lorsque Æ, décroît). 

La propagation de ces chocs pourra être ralentie par 
l’agitation thermique qui perturbe la rangée et par un 
désordre de masse. Ces effets sont faibles (= 1074 pour 
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Fic. 2. — Énergie perdue par collision d’après Silsbee [8] 


le désordre de masse dû au mélange isotopique du 
cuivre, mais cet effet peut être important pour les 
alliages). L'énergie sera surtout perdue par interaction 
avec les rangées voisines (voir fig. 2 pour le cuivre) et 
ceci limite la propagation du phénomène à une centaine 
de distances interatomiques. 


CONCLUSIONS DE SILBSEE. — 1; On a vu que, pour 
le cuivre, lorsque l’énergie a une valeur inférieure 
à 50 eV, l'effet de « canalisation » devient important ; 
or un tel processus ne donne naissance à aucun « dépla- 
cement » (c’est-à-dire interstitiel) supplémentaire. On 
peut alors estimer que la valeur de 50 eV représente 
une assez bonne évaluation de Æ4, minimum d’énergie 
nécessaire à « déplacer » un atome de cuivre. Cette 
valeur n’est pas en excellent accord avec les résultats 
expérimentaux [10] (æ 22 eV). Cependant cet argu- 
ment est très qualitatif, car un meilleur accord serait 
obtenu en modifiant légèrement le type d'interaction ; 
d’autre part, ceci suppose qu’il n’y a pas transport de 
matière, en même temps que d’énergie, au cours de la 
succession des collisions. 
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2. Un atome pourra être éjecté en position inters- 
titielle si la chaîne est interrompue, comme à la ren- 
contre d’une dislocation. Eventuellement le dernier 
atome de la chaîne pourra être éjecté de la surface du 
cristal si son énergie est suflisante, ce qui explique 
les expériences de Wehner. 

3. Une vibration anormale d’un atome pourra être 
transportée le long des rangées denses en un temps 
relativement court et à une distance assez grande, 
ce qui bouscule la théorie des « thermal spikes ». 


Extension du problème. — Selon les anciennes théo- 
ries, l'étendue du domaine perturbé par un atome pri- 
maire est de l’ordre de son libre parcours moyen. Des 
expériences d’élasticité [11] à basses températures 
semblent indiquer un ordre de grandeur plus élevé 
( 100 distances interatomiques). Leibfried [12] re- 
prend alors les idées de Silsbee pour rendre compte de 
ce fait. 

Au moment de la collision l’atome projectile se 
trouve en P (fig. 3). Comme la quantité de mouvement 


ETC: 


est transférée selon le vecteur joignant les centres des 
deux sphères dures, le deuxième atome sera projeté 
suivant PB, tandis que le premier partira dans une 
direction perpendiculaire (masses égales). Deux types 
de collisions peuvent alors se produire : 

a) Les collisions « canalisées » ne peuvent se produire 
que si « < 2, c’est-à-dire si R>d/2, ce qui correspond 
à une énergie limite £,4 au-dessus de laquelle ce pro- 
cessus disparaît. On voit (fig. 4) l’angle limite 6, pour 


Fic. 4. 


ce type de collisions et l’angle 0, au-dessus duquel 
il n’y a plus du tout de collisions. Comme l’atome pro- 
jectile est laissé dans une position plus proche de son 
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propre site que du site suivant, tous les atomes peuvent 


revenir à leur position initiale. Il se produit essen- 
tiellement un transport d’énergie (fig. 5). 


101] 


[001] 
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F1G. 5. — Canalisation des chocs. 


b) Si x > 2, les angles 0, vont aller en croissant, 
c’est-à-dire que les atomes vont s’écarter de plus en 
plus de la rangée ; c’est ce que nous appellerons le 
phénomène de « dispersion des chocs ». D’autre part, 
l’atome projectile est plus proche du site suivant, après 
le choc ; il est donc probable, dans une telle suite de 
collisions, qu’uné lacune se produise au début de la 
chaîne et qu’un interstitiel soit transporté le long de la 
rangée à une distance d’autant plus grande que l’angle 
initial sera petit. Leibfried appelle ce processus des 
« crowdion collisions » par analogie avec la vieille 
idée qui envisageait la’ possibilité d’interstitiel (statique) 
dans une rangée compacte, provoquant ainsi un « tas- 


[ooi] 


F1G. 6, — Dispersion des chocs. 


sement » des atomes voisins ; ici il s’agit évidemment 
d’un « crowdion dynamique » (fig. 6). 

Pour avoir des effets à longue distance, il faut que 
les quantités de mouvement fassent de petits angles 
avec la direction de la rangée. Leibfried poursuit alors 
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ses calculs en posant &« = 2 +1 avec [nl & 1 ce qui 
signifie que l’énergie du premier atome de la chaîne 
doit être de l’ordre de Ex. 

En développant les formules (1) pour de petits 
angles, on trouve que l’énergie du n-ième atome est : 


RE 2 
E, = E, F + 05 : 05 —| . 


On peut ainsi déterminer le nom bre maximum n*(1) 
de collisions dans un « crowdion » (en posant Æ, = 0 
et n > 0) et l’énergie maximum que peut transporter 
un « focuson » (7 co et n < 0). En fait le nombre de 
collisions est limité par l’énergie perdue dans chaque 
choc, principalement à cause de l’action des atomes 
voisins. Pour un potentiel du type de Huntington 
(avec des constantes choisies pour que £a æ 22 eV) 
et pour n = 0, Leibfried a trouvé que n*(0) æ 100. 

En se plaçant maintenant dans un cristal cubique à 
faces centrées, on peut évaluer, par des considérations 
sur les angles solides, la probabilité pour qu’un atome 
d'énergie £, = Ey donne naissance à un phénomène de 
«canalisation des chocs » ou de « dispersion des chocs ». 
Le calcul du nombre de collisions par atome primaire 
d'énergie Æ nécessite la connaissance du nombre de 
collisions S(E, Æ;) dE, dans lesquelles une énergie 
comprise dans l'intervalle (£,, Eo + dE5) est trans- 
férée à un atome secondaire. Utilisant un modèle dû 
à Kinchin et Pease, Leibfried trouve que le nombre de 
« crowdions » se propageant plus loin que n*(rn) est 
donné par : 


où y est le nombre d’interstitiels prévu par l’ancienne 
théorie. On voit que le nombre d’interstitiels pouvant 
être créé à grande distance est faible par rapport au 
nombre total d’interstitiels. 

Le nombre de « focusons » est au contraire beaucoup 
plus grand (20 °/, de v pour n* æ 50). Notons que, 
selon Leibfried, les collisions « canalisées » ne trans- 
portent que de l’énergie ; il n’y aura production d’in- 
terstitiels que lorsque la chaîne sera interrompue, par 
exemple près d’une dislocation, ce qui provoque son 
bloquage et fournit une explication des expériences 
d’élasticité précédemment citées. Selon ce modèle, 
Gonser [13] explique la dilatation de l’uranium-v sous 
irradiation. ù 


Résultats expérimentaux, — Thompson et Nelson 
[14, 15] ont bombardé des feuilles minces d’or, pré- 
sentant de fortes orientations préférentielles, avec des 
protons d’énergie 0,3 MeV en opérant par transmis- 
sion. Ils ont observé une éjection ,préférentielle d'ions 
selon les directions (110), la masse déposée sur le collec- 
teur'étant indépendante de l’angle d’incidence.Connais- 
sant le nombre d’atomes éjectés par proton, et utili- 
sant les méthodes de calcul de Leibfried ils trouvent 
que si Eno æ 800 eV,le parcours moyen d’un « focu- 
son » est de l’ordre de 350 À 
(d = 2,88 À), Bd — 15 et À — 0,8 x 105 eV 


(V(r) = À e—#) 


et que l’énergie perdue par choc est de l’ordre de 7 eV. 
Le potentiel de répulsion pour l’or est alors 


V(r) = 8 X 105 exp — 15 r/d 
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.ce qui conduit, en négligeant l’effet des électrons libres 


à une valeur de la compressibilité 0,33 X 10712 cm? 
dyne—1 en accord raisonnable avec la valeur expéri- 
mentale : 0,54 X 10—12 cm? dyne—, 

Par bombardement d’ions (A+, Xe+ énergie 10 keV), 
sur des couches d’or, de cuivre, et d’argent, ils ont 
observé un dépôt d’atomes neutres éjectés de la surface 
cristalline du même côté que le faisceau incident, selon 
trois directions préférentielles (110), (100) et (111). 

Pour la direction (110), ils interprètent le phéno- 
mène par une « canalisation » de l’énergie suivant la 
rangée. Suivant l’argument qualitatif de Leibfried, ce 
phénomène a lieu sans transfert de matière. 

Pour la rangée (100) il est clair (fig. 7) que les pre- 


miers voisins (atomes B) ont une influence prépon- 
dérante devant celle de l’atome suivant de la rangée. 
L’approximation des sphères dures n’est plus valable 
(paramètre. d’impact non négligeable) pour ces colli- 
sions à incidence rasante. On peut définir une énergie 
de focalisation £4 au-dessus de laquelle la déflexion ne 
permet pas à l’atome A, de rencontrer l’atome A, 
et une énergie minimum Æ;, au-dessous de laquelle 
l’atome À, ne peut traverser l’anneau des atomes B. 
Le même phénomène a lieu pour la rangée (111). La 
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canalisation des chocs a lieu ici avec transport de ma- 
tière. 


Er (eV) E; (eV) 
DIRECTIONS. Cu Au Cu Au 
(110) 60 800 1 1 
(100) 65 700 10 90 
(111) 490 7 300 AA 500 


Valeurs calculées avec 
V(r) = 2 X 104 X exp (— 13 r/d) pour Cu [12] 
et V(r) = 8 X 105 X exp (— 15 r{d) pour Au. 


Ces résultats sont à comparer à ceux de Koedam [16]. 
En bombardant avec des ions (A+. Kr+, Ne+, 50 à 
2 000 eV) sous incidence normale aux plans (100) des 
cristaux de cuivre, il a observé 4 directions (110) d’éjec- 
tion préférentielle ainsi que la direction. (100) normale 
à la surface bombardée ; l'intensité relative de ce 
dernier dépôt croît avec l’énergie des ions incidents, 
ce qui semble bien indiquer la présence d’un seuil mi- 
nimum Æ1%, 


a 


Caleuls numériques de Gibson, Goland, Milgram et 
Vineyard [17]. — Considérant que les méthodes de 
calcul analytique ne sont pas adaptées au problème 
si on veut tenir compte des interactions entre des 
atomes répartis régulièrement dans un espace à trois 
dimensions, cette équipe a lancé dès 1957 une série de 
calcuis numériques, grâce à l’utilisation d’une machine 
IBM 704. 

Les calculs ont été faits en prenant pour modèle de 
petits cristallites de cuivre (de 500 à 1 000 atomes) 
supposés isolés à l’intérieur d’un grand cristal, à la 
température de 0 °K pour ne pas avoir à tenir compte 
de l’agitation thermique. Les forces utilisées sont de 
trois types : 

a) L’interaction entre les noyaux a la forme clas- 
sique de Born-Mayer : 

V(r) = 4 exp [- 8] 


la plupart des calculs a été faits avec : 
A=0,051,  p = 13. 


Des essais effectués avec d’autres couples de valeurs 
n’ont amené aucune différence appréciable dans les 
résultats statiques et dynamiques des calculs. La seule 
différence se trouve dans l’estimation de Æ4. Le poten- 
tiel choisi conduit à £4 25 eV dans la direction (110) 
de 25 à 30 eV dans la direction (100) et æ 85 eV dans 
la direction (111). 

b) Les forces de cohésion (dues aux électrons de con- 
duction) sont représentées par une force constante 
appliquée aux atomes de la surface et dirigée vers 
intérieur du cristal. 

c) L’action des atomes de la matrice cristalline sur le 
cristallite considéré est représentée par une force élas- 
tique proportionnelle aux déplacements des atomes et 
une force dissipative proportionnelle à la vitesse des 
atomes de surface. En faisant varier la position de 
l’atome mis en mouvement et la taille du cristallite, on 
n’a pas observé de différences notables dans les dom- 
mages résultants, ce qui justifie le choix de ces forces. 
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Un atome quelconque de la maille est mis en mouve-. 
ment avec des énergies Æ, variant de 20 à 400 eV, 
dans des directions faisant de petits angles avec les 
rangées compactes et une succession de « passages » 
dans la machine permet de déterminer la position de 
tous les atomes en fonction du temps. 

— Le phénomène de «canalisation des chocs » a lieu 
pour des énergies £, < 30 eV dans la direction (110) et 
< 40 eV dans la direction (100). Contrairement aux 
hypothèse de Leibfried on peut avoir, dans ce modèle 
un transport d’interstitiel aussi bien que l’énergie si 
l'énergie initiale £, est supérieure à Æa. Pour des éner- 
gies supérieures aux valeurs limites Er, On voit appa- 
raître le phénomène de « dispersion des chocs ». Les 
défauts créés sont des paires de Frenkel avec des 
lacunes situées près du point de départ des chaînes de 
collisions et des interstitiels dédoublés (split confi- 
guration), seuls stables. 

—— Concentrons-nous plus spécialement sur la 
rangée (110). On peut définir un paramètre de canali- 
sation : 

A(Eï) = (Oi+3) JO. 


La série de calculs faits montre que À dépend très 
peu de 6 et l’on peut tracer une courbe moyenne A(F;) 
(fig. 8) ; d’autre part on peut écrire que l’énergie perdue 


0 50 100 150 
ENERGIE (eV) 4 


Fic. 8. 


par choc dans la rangée (110) obéit approximativement 
à la relation : 


AE; = [(2/3) + Eo sin? (Or + Oi421)] eV (2) 


formule valable pour de faibles angles et pour les 
domaines d'énergie considérés. 

Si la chaîne se disperse, 0 croît et l’énergie trans- 
portée décroit plus rapidement jusqu’à ce qu’elle re- 
tombe au voisinage de Æ4. La dispersion des chocs 
entraîne évidemment la mise en mouvement des ran- 
gées voisines. 

Si £, > Ea, nous aurons un transport d’interstitiel. 
Les calculs montrent qu’un interstitiel est créé dès que 
l’énergie transportée tombe à 3,5 eV. La connaissance 
de À et la relation (2) permettent alors de prévoir le 
nombre de remplacements produits dans la rangée et 
par conséquent le nombre de distances interatomiques 
que peut parcourir un « focuson » (fig. 9). 

La figure 10 illustre une énergie de calculs numé- 
riques. On communique au temps t — 0, une énergie 
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Fic. 9. — Les abaques représentées ici permettent de 
déterminer en fonction de 6, et de Æ, le nombre de dis- 
tantes interatomiques parcourues par un interstitiel 
« canalisé » le long de la rangée (110) [17]. - 

Les courbes s’arrêtent pour E © 30 eV. Pour des 
énergies inférieures nous n’avons plus en effet que de 
« l’énergie canalisée » et le calcul ne s’applique pas à ce 
cas. 


de 40 eV à l’atome A ; sa vitesse initiale fait un angle 

de 220 1/2 avec l’axe y (100) tout en restant dans le 

plan représenté. On voit nettement un phénomène de 

canalisation des chocs avec transport d’interstitiel le 

long de la rangée AC (100). Dans la direction AB, au 

contraire, il y a seulement canalisation de l’énergie sauf 
eut-être l’atome suivant À qui semble vouloir pro- 
uire un interstitiel dédoublé. 

On peut hésiter dans le choix des hypothèses de 
départ mais le calcul numérique offre de grandes possi- 
bilités. Un des aspects attachants du travail de 
Vineyard et son équipe est de permettre de suivre 
i’évolution du phénomène au cours du temps ; un film 
illustre d’ailleurs leurs résultats et donné une vision 
d'ensemble des mouvements simultanés des atomes de 
la maille. 
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F1G. 10.— La figure représente un plan (100). Les positions 
d’équilibre des atomes du plan sont entourées d’un 
cercle. Les gros points noirs indiquent la position des 
atomes au temps t — 324.107 1% secondes. Les petits 
points représentent Les positions d’équilibre des atomes 
des plans adjacents, les croix leur position finale. 
1 atome À a une énergie initiale de 40 eV [17]. 


Conclusion. — Les travaux récents sur la propa- 
gation des chocs dans les solides permettent une meil- 
leure compréhension des dommages produits par irra- 
diation. Par exemple, ils montrent la possibilité d’une 
large séparation lacune-interstitiel et ils apportent des 
lumières sur le mécanisme du désordre dans les alliages. 
Malgré certaines divergences, les modèles proposés 
conduisent à une même description qualitative du phé- 
nomène. De plus ils ont tenu compte de l’influence du 
réseau sur le comportement des atomes et l’on peut 
espérer que la prise en considération de ce fait fonda- 
mental amorce une théorie plus cohérente des défauts 
créés par irradiation dans les solides. 


Manuscrit reçu le 13 décembre 1960. 
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LETTRES A LA RÉDACTION 


RÉACTIONS (p,p an) SUR LE CÉSIUM 


par CG. Rien 
Laboratoire de Physique Nucléaire, Orsay (Seine-et-Oise). 


Dans le cadre de l’étude des réactions de spallation 
actuellement en cours au Laboratoire de Physique 
Nucléaire d'Orsay [1],, nous avons irradié du césium 
par le faisceau de protons du synchrocyclotron pour 
des énergies variant entre 30 et 155 MeV et nous avons 
étudié les sections efficaces de formation de divers 
isotopes du césium [2]. 

Ces divers isotopes se désintègrent essentiellement 
par capture électronique sur la couche X. Nous avons 
observé les isotopes suivants : 


Y CARACTÉRISTIQUES 


ISOToPE PÉRIODE 
B2CS (LS 670 keV 
ICS SLT EAU — 
CS 31e 380 keV 
(ES 62h 410 keV 


La détermination des sections efficaces se fait par 
mesure de l'intensité du rayonnement XX de 30 keV 
émis avec décomposition de la courbe de décroissance 


133 132 


en mb Cs'{(. p pn à €s 


x s théorique 
le 
i T exp 


100 E 


Error 


d’une part et par mesure de l'intensité des rayonne- 
ments y caractéristiques d’autre part. 
Le césium est irradié sous forme de poudre de chlo- 


rure de césium. À chaque cible, on associe un moniteur 
de graphite ou d'aluminium suivant l'énergie. 

Le césium est ensuite purifié par séparation chi- 
mique et les mesures d'intensité des divers rayonne- 
ments sont effectuées par des méthodes de spectro- 
graphie y en utilisant des cristaux de INa comme 
scintillateur. Les impulsions sont analysées par un 
sélecteur à 60 canaux. 


Ce 133 


AS: 
Cs { p,p2n) Cs 
le 
} S exp- 


AC: thæorique 


50 100 


HTG:22. 


Nous avons tenu compte de diverses corrections 
dans le calcul des sections efficaces : phénomènes 
d’autoabsorption pour le rayonnement de 30 keV, 
rapports capture X/capture Z, coefficients de conver- 
sion interne et influence des filiations radioactives, 
ce qui nous à permis de déterminer la section eflicace 
de la réaction (p, 5n) en même temps que celle de la 
réaction (p, p 4n), mais, d’autre part, de n’observer que 
la section efficace cumulative (p, p 6n) + (p, 7n). 

Les erreurs ont été évaluées à 15 % pour la réaction 
(p, pn), 25 % pour les réactions (p, p 2n) et (p, p 6n) 
et 50 % pour les réactions (p, p 4n) et (p, 5n) (voir 
figures 1 à 4). , 

S1 on se réfère au schéma classique des réactions de 
spallation, processus en deux étapes :interaction directe 
et évaporation (en négligeant l’évaporation des proz 
tons), on peut interpréter l’allure des courbes obtenues. 
La position à 60 MeV du maximum de la réaction 
(p, pa) correspond à 45 MeV'environ d’énergie de eut- 
Of pour le proton, augmentés d’environ 15 MeV, 
énergie à laquelle un neutron a la plus grande proba- 
bilité d’être évaporé. On observe ensuite, pou: les 
autres réactions, un déplacement du maximum d’envi- 
ron 10 MeV par neutron supplémentaire évaporé. 
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D'autre part, la valeur notable des sections efficaces 
à plus haute énergie peut être attribuée à l'interaction 
directe. 


DH  stppanl + 07 &ip,5n) 
ä ® 


en mb 


ip,p4n) +1,09 {(p,5n} 
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| Fic. 3. 
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Nous avons effectué des calculs théoriques à partir 
des résultats de Métropolis [3] pour l'interaction di- 
recte, et Jackson [4] pour l’évaporation. De façon 
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générale, nous constatons une assez bonne concor- 
dance à 155 MeV entre théorie et expérience, mais les 
calculs théoriques ne rendent pas bien compte des 
maxima des différentes courbes, sauf pour la réaction 
(p, pn) pour laquelle nous observons un écart d’un 
facteur 4 environ. Ce fait est commun aux réaction 
(p, pa) et laisse prévoir que les approximations effec- 
tuées ne sont peut-être pas toutes justifiées et qu'il 
peut exister un autre processus dont on ne tient pas 
compte. 
Lettre reçue le 22 juin 1961. 
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SUR LE MÉCANISME 
DE LA RÉACTION 12C(y, 3x) à 17,6 MeV 


Par Wilfrid Sésaoux et Henri GAUVIN, 


Institut du Radium, Laboratoire Joliot-Curie 
de Physique Nucléaire, Orsay. 


Nous avons utilisé pour cette étude des émulsions 
ionographiques diluées en gélatine et le rayonnement 
y de 17,6 MeV de la réaction ‘Li(p, y). Dans un précé- 


N/0,5 MeV 


— 200 


| L 400 


(e] D 10 
F1G.°14 — "Distribution de F*, 


dent travail [1], sur la section efficace de la réaction 
2((y, 3x), sont données des précisions sur les condi- 
tions de l’expérience et les mesures. 

Le mécanisme de la réaction 12C(y, 3x) a été discuté 
antérieurement par plusieurs auteurs [2-10]. Les 
processus envisagés sont 1° la tripartition directe du 
noyau de 2?C, et 20 l'émission par le noyau excité de 
2C d’une particule &,, suivie de la désintégration du 
noyau résiduel de $Be en deux particules &, et æs. Un 
moyen de vérifier cette dernière hypothèse consiste à 
examiner la distribution des énergies d’excitation du 
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8Be. À chaque trace observée, supposée produite par la 
particule «, correspond une énergie d’excitation 
possible E*,. Les niveaux du 8Be excités doivent faire 
apparaître des pics dans la distribution de E*. Évidem- 
ment, comme on ne peut identifier «&,, un tiers seule- 
ment des valeurs de E* calculées donnent les niveaux 
excités du $Be. La figure 1 montre la distribution 
obtenue. Cette distribution présente un pic corres- 

ondant au niveau connu de 2,9 MeV, de largeur 
1,2 + 0,3 MeV [11], et un pic vers 8,7 MeV. Elle est 
compatible avec l’interprètation donnée de leur dis- 
tribution par Glätthi et al. [8], qui ont admis, outre 
l'intervention du niveau de 2,9 MeV, celle de deux 
niveaux à 4 MeV et 9 MeV, et d’un niveau éventuel 
vers 7,5 MeV. 

. Cependant, en dehors du fait que l’existence de ces 
niveaux autres que celui de 2,9 MeV est incertaine [11], 
il faut remarquer que cette analyse est insuffisante. En 
effet, on a, en négligeant la quantité de mouvement du 
photon : 


E*=E,—Q—{1+ m|M) Ex, 


où Æ, — 17,6 MeV est l'énergie du photon, Q’ le seuil 
de la réaction ?C(y, «)$Be, m la masse de 4He et M 
la masse de $Be. 

Les pics de la distribution de Æ* reproduisent donc 
ceux de la distribution de Æ, (fig. 3). Telegdi [3] a 
montré qu’on peut rendre compte de la distribution 
de E, observée si l’on suppose que la réaction 1C(Y, 3x) 
à 17,6 MeV implique seulement le niveau de 2,9 MeV 
du Be, à condition de tenir compte de la corrélation 
angulaire entre «&, et 2,4, calculée en admettant que 
l'absorption du rayonnement y est une transition 
mixte 0,5 M1 + 0,5 Æ2, avec peut-être une faible pro- 
portion de Ef. 

Nous avons repris la méthode proposée par Telegdi 
pour analyser nos résultats. 

Le niveau de 2,99 MeV du fBe est un niveau 2+, 
T = 0[11]. Pour une transition y du type E2 on a la 
séquence : 


12C(j = 0+) + hv —> 12CX(j = 2+) 
—+ 8Be*(j = 2+) + où (I). 


On a principement ! — 0 et il n’y a pas de corréla- 
tion. 
Pour une transition M1 on a la séquence : 


1 (j = 0) + y => 1C#(j = 1+) 
> °Be*(y.= 2F) LG 2), 


Le niveau de 2,9 MeV de Be se désintègre par 
émission d’une onde de ! — 2. La fonction de corréla- 
tion entre «, et &,,:, donnée par Telegdi est : 


W(0) = (15/4) cos?0 (1 — cos? 0). 
Pour une transition £1 on a la séquence : 
PC( = 0+) + hv — PCA (7 = 1—) 
ne 2 Nu Bt AN 


On a principalement ! — 1 et la fonction de corréla- 
tion correspondante, calculée à l’aide des formules 
de Lloyd [12], est 


W(8) = (1/2) (1 + 3 cos? 0). 
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Les distributions de E,; et Euo,s correspondant aux. 


différents cas ont été calculées en tenant compte de la 
largeur L — 1,2 MeV du niveau de 2,9 MeV du fBe. 
Elles sont représentées sur la figure 2. 


p (E) 


absorption M1 


Ex 


absorption E 2 


Ex 


absorption E1 


Fic. 2. — Distributions théoriques de Æ,. 


Dans le cas d’une tripartition directe de 22C*, en l’ab- 
sence de corrélation entre les particules émises, la = 


distribution de Æ, est donnée par : 


p(Es) = cte X EX? [Er — (3/2) EJ4?, 


où Æy est l'énergie totale des particules « [5]. Cette - 


relation est représentée par une demi-ellipse. 


La figure 3 montre la distribution expérimentale - 


de F, et deux distributions calculées. 

La première courbe correspond à l’hypothèse de 
Telegdi d’une absorption du rayonnement y M1 + Æ2, 
Le meilleur accord avec la distribution expérimentale 
est obtenu pour le mélange 0,6 M1 + 0,4 E2. La 
2€ courbe correspond au cas d’une transition E1 pure, 
C'est le cas qu’il faut envisager si l’on admet que le 
niveau large de 17,22 MeV de 2C(T — 1,27 MeV, 
1 —1—, T = 1) est responsable de l'absorption du 
rayonnement y de 17,6 MeV [13]. 


L’accord avec la distribution expérimentale est : 


moins bon que pour la transition M1 + Æ2, mais un 
déficit d'étoiles comportant une particule &« d'énergie 


< 0,5 MeV (de parcours < 1,5 1) pourrait peut-être . 


s’expliquer par le fait que de telles étoiles sont très 
difficilement identifiables. 

. Toutefois, cette hypothèse rencontre une difficulté 
signalée par Goward et Wilkins [9] : 


Non 


r 


mu INo 11 


le NS LAC dR ES 
Es 


La réaction 2C*(T = 1) — 8Be*(T = 0) + &(T = 0) 
n’est possible qu’à condition d'admettre une certaine 
ue des états de spin isotopique de 2C* ou 

e*, 


N/0,25 MeV 


Spectre observé 
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F1c. 3. — Distribution expérimentale 
. et distributions calculées de £E,. 


Un moyen de préciser les conclusions de cette ana- 


. lyse serait d’étudier directement la corrélation angu- 


laire entre &, et «4. Malheureusement, la dispersion 
de Æ,, due à la largeur du niveau de 2,9 MeV du 
8Be rend très hasardeuse l'identification de la parti- 
cule œ. 

Nous remercions Mme Boudét et Me Dormeau pour 
l’aide qu’elles nous ont apportée. 
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QUELQUES PROPRIÉTÉS 
DE GRENATS COMPLEXES DE TERRES RARES 


par Gérard VizLers et Jean Lorters, 
Laboratoire de Magnétisme et de Physique du Solide 


et Laboratoire des Terres Rares, C. N. R. S$. 
Bellevue (Seine-et-Oise). 


Dans les travaux antérieurs [1, 2, 3] nous avons 
exposé les propriétés de grenats mixtes de terres rares 
répondant à la formule : 5 Fe,0,. æ Yt,0:, (3 — x) 
Yt,0, (Ytet Yt' désignant soitl’yttrium, soit une terre 


42,450 


42,400 


Fret 


rare. Nous donnerons ici les résultats expérimentaux 
obtenus sur des grenats plus complexes représentés 
par la relation moléculaire : 5 Fe,O,, x Yt,0:, 
2Yt' Os, [8 — (x + 2)], Yt',0, (Yt, Yt!, Yt” désignant 
soit l’yttrium, soit des terres yttriques). A cet effet 
nous avons choisi deux groupes distincts : le premier 
comprend l’yttrium et deux terres yttriques, le second 
comporte trois terres yttriques différentes dans la 
maille du grenat. Comme exemples, caractérisant ces 
deux groupes, nous avons retenu les systèmes : 


5Fe203, ZY 303, 2Gd203, [3 — (x + z)] D y°03 et 
5Fe:O3, xGd203, 2DY203, [8 — (x + z)] Er:O», 


dans lesquels nous avons fait varier x et z de 0 à 3 
avec (æ + 2) < 3 et en nous limitant à la condition 
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TABLEAU 1 
5Fe:03, 2Y203, 2Gd203, [3 — (x + 2)] Dy203. 
0c0K 0 °K 
GATAO zGd203 [3 TT (& + z)] Dy203 dy À Ld To À Of 0K Pre A 
2.0 0,5 0,5 12,395 1,0% 555 — 51 
4,0 1,0 4,0 12,412 1,08 557 155 152 
0,5 0,5 2,0 12,402 1,07; 558 193 191 
0,5 2,0 0,5 12,435 1,09; 560 242 212 
TABLEAU 2 
5Fe:03, 2Gd203, 2Dÿ203, [8 — (x + z)] Er2O3. 

: OSROIE OK: 
azGd,03 zDy:0s [3 — (x + z)] Er,O; & À ro À 0 °K NES Éd ME 
2.0 0,5 0,5 12,434 1,09, 599 234 236 
4,0 1,0 1,0 12,403 1,0% 596 190 191 
0,5 0,5 2,0 42,372 4,05, 550 130 125 
0,5 2,0 0,5 12,400 1,0% 55% . 208 203 


2z — 3 — x. On remarquera d’une part que si x et z 
s’annulent en même temps nous sommes ramenés à 
l’étude de grenats simples, d’autre part que si æ — 0 
(ou bien z — 0) nous sommes dans le cas des systèmes 
de grenats mixtes binaires. Des points expérimentaux 
obtenus antérieurement [1, 2], concernant ces deux 
conditions particulières, ont été reportés sur les 
courbes des figures 1 et 2. 


OPK 
[ 
À 
5e0 À 
?— 
Pt 
540 | 
18,300 1,350 42,00 12,€0 A À 
ICRA 


Ces différentes compositions ont été préparées par 
coprécipitation à partir d’oxydes de terres rares et de 
fer de pureté supérieure à 99,9 % en utilisant les 
méthodes déjà décrites [1]; les propriétés physiques 
mentionnées sont relatives à des matériaux frittés à 
1350 0C pendant quatre heures dans l’air et refroidis 
en quinze heures environ dans le four. 

Sur ces ferrites nous avons mesuré les valeurs des 


paramètres de la maille cristalline à,, des températures 
de Curie 0, °K et de compensation 60. °K ; nos résultats 
sont consignés, pour chacun de ces deux groupes, dans 
les tableaux { et 2. 

La courbe de la figure 1 représente la variation de 
a, en fonction du rayon ionique moyen 7, des ions 


terres rares entrant dans le réseau du grenat (r, à été … 


calculé à partir des dimensions d’ions données par 
Goldschmidt et de leurs proportions relatives). Il 


ressort de cette courbe que le paramètre croît régu- : 


i 
$ 
# 
#5 
j 
Li 
! 
DT 
: 


lièrement et dans le même sens que r,. Sur la figure 2 


nous avons reproduit la courbe montrant la variation 
des points de Curie 0; 0K en fonction de a,. On re- 
marque, comme nous l’avions déjà constaté sur des 
grenats mixtes binaires [1, 2], que les points de Curie 
sont d'autant plus élevés que le paramètre cristallin 
est plus grand. Dans les tableaux 1 et 2, pour les com- 
positions étudiées, nous mentionnons également les 
températures des points de compensation 60,°K théo- 
riques. Ces dernières ont été obtenues graphiquement 
à partir des aimantations spontanées des trois sous- 
réseaux 164, 24d et 24c calculées par R. Pauthenet [4], 
en fonction de la température ; nous avons tenu compte 
évidemment de la proportion relative des différents 
ions terres rares placés sur les sites 24c et nous avons 
adopté le modèle ferrimagnétique proposé par 
L. Néel [5]. On notera que les résultats théoriques et 
expérimentaux sont en bon accord. 


Lettre reçue le 21 juillet 1961. 
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UNIFORMISATION 
DU VECTEUR CHAMP MAGNÉTIQUE 
SUR L’AXE D’UNE HÉLICE 


Par M. PaPpouLaRr, 


Laboratoire d'Électronique et de Radioélectricité, 
Fontenay-aux-Roses (Seine). 


1. — Il est intéressant, en métrologie par exemple, 
de réduire le champ magnétique sur l'axe d’un solé- 
noïde (hélice) à sa composante longitudinale, qui seule 
est de répartition axiale uniforme (les méthodes 
modernes de mesure par résonance magnétique nu- 
cléaire portent sur le vecteur H total). Nous supposons 
qu'il s’agit d’un solénoïde infiniment long (autrement 
dit : les considérations qui suivent ne nous dispensent 
pas des corrections de bout). 

2. — Le calcul permet de déterminer de façon rigou- 


_reuse la structure de ce champ [1]: a) le vecteur 


champ, en un point M de l’axe a une longueur cons- 
tante et fait avec l’axe un angle constant ; b) ce vec- 
teur est situé dans un plan perpendiculaire à MA, 
À étant le point du solénoïde, qui a la même cote 
que M. 

Compte tenu de ces résultats, on peut améliorer 


 luniformité du champ au voisinage de l’axe, en inter- 


calant dans la première hélice, une deuxième, iden- 


tique, parcourue dans le même sens par le même cou- 


rant, mais décalée d’un demi-pas par rapport à la 
première [2]. 

Sur l'axe, on obtient ainsi une compensation à la 
fois simple et rigoureuse des composantes transversales. 


En effet, À, et A, étant symétriques par rapport 
à Oz, H,, H, et Oz sont dans un même plan, et H, 
étant équipollent à H;, ne peut qu'être symétrique 
de H, par rapport à Oz. Il est intéressant de noter que 
la compensation sur l’axe est rigoureuse, quel que soit 
le pas commun des deux solénoïdes associés. Dans Ja 
pratique,on pourra enrouler en hélice un cordon bifilaire 
dont on mettra les deux brins en série. On se heurte 
évidemment à des difficultés technologiques quand on 
veut atteindre la haute précision requise en métrologie, 
notamment sur la longueur du pas [3]. 
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3, — Pour les applications, il est nécessaire d’éva- 
luer le degré de compensation obtenu aux environs de 
l'axe, afin de délimiter le cylindre utile où le taux d’uni- 
formité désiré est réalisé. Reprenons le calcul du champ, 
mais cette fois en un point quelconque M(X, Y, 2). 
Nous ferons au voisinage de l’axe, un développement 
en série de l’expression trouvée, que nous limiterons 
au 1er ordre. 

Ainsi par exemple pour la composante 72,, l'élément 
de cote z — Àt de la {re hélice, donne : 


3R(X cos t + Y sint) 
R? + (Z — ki) 


el 
R(Z — ki) cos t + Rksint 


(dHaha = (4H) [1 + 


L'élément de même cote de la 22 héiice, donne : 


3R(X cos t + Ysint) 
R? + (Z — kt)? 


(Has = (A), [1 


kY | 


. R(Z — kt) cos t + Rksint 


(d;), représente la contribution de l'élément z — ki 
de la 17e hélice au point M(O, O, Z) sur l'axe : 


.R(Z — kit) cos t + Rksint 


(dy) = 1 RE ee MÈEE di. 


On opère le changement variable : 
Z — kt = — Ru. 


Et on obtient au total, après des calculs simples : 


His le nl | 
ee er 


77 OO 


k R 
Gt 07 re = cos 2 ee 
+R sn k R k : 
À Res 


Et Ru k Ru 
TT - DR LES 
u sin k + R °° see 
â (1 u2)5/2 


Il est commode de poser : « — 2R/k. 
On utilise alors la relation : 


ee cos au 27 œ\® 
7e K 
fe RAT) CEE USE ET ONE TRE (:) n(o) 


où X, est la fonction de Bessel modifiée, de 2€ espèce. 
On obtient : : 


Gi A 
Dire Y cos 


Ce 


Fa Le im (? Æ{a)) — = im (x 0) 
21 27 27 . 
1 [y sin 22 y cos 22 te? X D RE 
ek (4 sin k Y cos A [a? Kilx) + 20 Kolæ)] 


Le comportement des fonctions X,(&) quand « tend 
vers Zéro, fait que le premier terme entre crochets est 
nul. Quant au second, on peut profiter de la relation 
de récurrence : 


Ey+a(o) = Ayo) + ŸY Ava) 
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pour le mettre sous la forme : 
20 Hola) + (4 + «?) Hilo). 


Le calcul de 7, est mené suivant les mêmes consi- 
dérations. Nous négligeons la variation de 77, au voisi- 
nage de l’axe, qui est nulle au 1° ordre [4]. On arrive 
finalement à : 


: EN ZEN 2Z 

He HTC niG{o) sn — Aus | 
2Z F2 
H, = — LT ni G{o) É Poe + RE | 


avec : 


Ga) = 2x Kola) + (4 + à?) Kila). 


On voit immédiatement que lorsqu'on déplace le 
point M parallèlement à l’axe, à la distance r de celui- 
ci, l’extrémité du vecteur champ décrit une hélice de 
pas 2x k/2 et de rayon : 


kr nG(a)r[R. 


En l’absence de compensation au contraire, l’extré- 
mité du vecteur champ, en un point M, de l’axe, décri- 
rait une hélice de même pas 2x4 que l’hélice solénoïde 
et de rayon [1]: 


semi aves #9 = 2 (95) 


= Vo? + y? 


PrG:"2, 


Le degré de compensation réalisé au voisinage de 
l'axe, sera donc mesuré par : 
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AH transv. r 
H transv.  R F(«f2) 


Pour & _ 15, on a l'expression asymptotique : 


AH transv. 
H trans v. 


ei ne. 


« théoriques » dont les spires seraient sans ae d 
Il donne cependant l’ordre de grandeur de l’homo=" 
généité de champ réalisable dans un solénoïde (double) À 
de type ordinaire. 
Par exemple pour un solénoïde du genre du 
NBS 590 80 : 
2h —=2801mm, p — 2746 — Tim, 07m 0 3 
« étant voisin de 2 000, on voit que l’homogénéité 
est pratiquement parfaite. Par suite de la compen- 
sation des composantes transversales, le solénoïden 
double peut se construire avec un diamètre plus petit, 
sans nuire à l’uniformité du champ : on gagne alors 
sur l’importance des corrections de bout à effectuer. ki 
En application des résultats précédents, on peut 
envisager une méthode de compensation « naturelle » 
des composantes transversales du champ magnétique. 
Il suffirait d’enrouler en hélice une bande métallique, des. 
façon que l'intervalle entre spires soit faible devant law 
largeur de la bande. Pratiquement en effet, à chaque“ 
hélice-ligne de courant est alors associée une hélice 
«compensatrice » décalée d’un demi-pas. Il s’agit d’un un 
raisonnement analogue à celui que l’on fait en optique, a 
en ce qui concerne les fentes de diffraction. | 1] 
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[4] L'expression de Æ, comporte, de même que cales | 
de A, et H,, un terme indépendant de « : 


2 
PTE Mie [a Kyo) — ka? Æo(æ)]. 


De même que pour A, et H,, ce terme est nul, 
si bien que l’on peut effectivement négliger la varia= 
tion de /Z,. 


